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Abstract 

We show that eigenvalues of the family of Baxter Q-operators for supersymmetric integrable spin chains 
constructed with the gl(K\M )-invariant R- matrix obey the Hirota bilinear difference equation. The nested 
Bcthe ansatz for super spin chains, with any choice of simple root system, is then treated as a discrete 
dynamical system for zeros of polynomial solutions to the Hirota equation. Our basic tool is a chain of 
Backlund transformations for the Hirota equation connecting quantum transfer matrices. This approach 
also provides a systematic way to derive the complete set of generalized Baxter equations for super spin 
chains. 
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1 Introduction 



1.1 Motivation and background 

Supersymmetric extensions of quantum integrable spin chains were proposed long ago [1, 2] but the proper 
generalization of the standard methods such as algebraic Bethe ansatz and Baxter TQ-relations is still not 
so well understood, as compared to the case of integrable models with usual symmetry algebras, and still 
contains some elements of guesswork. 

Bethe ansatz equations for integrable models based on superalgebras are believed to be written according 
to the general "empirical" rules [3] applied to graded Lie algebras. Accepting this as a departure point, one 
can try to reconstruct other common ingredients of the theory of quantum integrable systems such as Baxter 
relations and fusion rules. For the super spin chains based on the rational or trigonometric i?-matrices, the 
algebraic Bethe ansatz works rather similarly to the case of usual "bosonic" spin chains, but the Bethe ansatz 
equations for a given model do not have a unique form and depend on the choice of the system of simple 
roots. The situation becomes even more complicated when one considers spins in higher representations 
of the superalgebra, and especially typical ones containing continuous Kac-Dynkin labels. The systematic 
description of all possible Bethe ansatz equations and TQ-relations becomes then a cumbersome task [4] -[6] 
not completely fulfilled in the literature. To our knowledge, a unified approach is still missing. 

In this article we propose a new approach to the supersymmetric spin chains based on the Hirota-type 
relations for quantum transfer matrices and Backlund transformations for them. Functional relations between 
commuting quantum transfer matrices are known to be a powerful tool for solving quantum integrable models. 
They are based on the fusion rules for various irreducible representations in the auxiliary space of the model. 
First examples were given in [7, 8]. Later, these functional relations were represented in the determinant 
form [9] and in the form of the Hirota bilinear difference equation [10, 11]. In [12]-[15] it was shown that 
transfer matrices in the supersymmetric case are subject to exactly the same functional equations as in the 
purely bosonic case. The Hirota form of the functional equations has been proved to be especially useful 
and meaningful [16, 17]. It is the starting point for our construction. 

The Hirota equation [18] is probably the most famous equation in the theory of classical integrable 
systems on the lattice. It provides a universal integrable discretization of various soliton equations and, at 
the same time, serves as a generating equation for their hierarchies. In this sense, it is a kind of a Master 
equation for the theory of solitons. It covers a great variety of integrable problems, classical and quantum. 

In our approach, quantization and discretization appear to be closely related in the sense that solutions 
to the quantum problems are given in terms of the discrete classical dynamics. What specifics the problem 
are the boundary and analytic conditions for the variables entering the Hirota equation. In applications to 
quantum spin chains, these variables are parameters of the representation of the symmetry algebra in the 
auxiliary space. For representations associated with rectangular Young diagrams they are height and length 
of the diagram denoted by a and s respectively. In the case of spin chains based on usual (bosonic) algebras 
gl(K), the boundary conditions are such that the non- vanishing transfer matrices live in a strip < a < K 
in the (a, s) plane [16, 17]. In the case of superalgebras of the type gl(K\M), the strip turns into a domain 
of the "fat hook" type presented in Fig. 1. 

In the nested Bethe ansatz scheme for gl(K), one successively lowers the rank of the algebra gl(K) — > 
gl(K— 1) (and thus the width of the strip) until the problem gets fully "undressed". This purely quantum 
technique has a remarkable "classical face" : it is equivalent to a chain of Backlund transformations for the 
Hirota equation [16, 17]. They stem from the discrete zero curvature representation and associated systems 
of auxiliary linear equations. This solves the discrete Hirota dynamics in terms of Bethe equations or the 
general system of Baxter's TQ-relations. 

The aim of this article is to extend this program to the models based on superalgebras gl(K\M). In this 
case, there exist two different types of Backlund transformations. One of them lowers K while the other one 
lowers M. The undressing goes until the fat hook is collapsed to K = M = 0. The result of this procedure 
does not depend on the order in which we perform these transformations but the form of the equations does. 
Different orders lead to different types of Bethe ansatz equations and Baxter's TQ-relations associated with 
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Figure 1: The domain ("fat hook") of non-vanishing transfer matrices T(a, s, u) for the super 'symmetric spin 
chain with the gl(K\M) symmetry. 

Cartan matrices for different systems of simple roots. In this way, the abundance of various Bethc equations 
and TQ-relations is easily explained and classified. All of them are constructed in our paper. Instead of 
K + 1 Baxter's Q-functions for the bosonic gl(K) algebra we recover (K + 1)(M +1) Q-functions (some of 
them arc initially fixed by the physical problem). More than that, we establish a new equation relating all 
these Q-functions which is again of the Hirota type. This "QQ-rclation" opens the most direct and easiest 
way to construct various systems of Bethe equations. Similar relations hold for the transfer matrices at each 
step of the undressing procedure. 

Our construction goes through when observables in the Hilbert space of the generalized spin chain are 
in arbitrary finite dimensional representations of the symmetry (super) algebra. We can also successfully 
incorporate the case of typical representations carrying the continuous Kac-Dynkin labels, as it is illustrated 
by examples of superalgebras gl and gl(2\l). 

Some standard facts and notation related to superalgebras and their representations are listed in Ap- 
pendix A. For details see [19]-[22]. Throughout the paper, we use the language of the algebraic Bethe ansatz 
and the quantum inverse scattering method on the lattice developed in [23] (see also reviews [1, 24] and 
book [25]). On the other hand, we employ standard methods of classical theory of solitons [26] and discrete 
integrable equations [27]-[31]. 

1.2 A sketch of the results 

We consider integrable generalized spin chains with gl(K |M)-invariant i?-matrix. The generating function 
of commuting integrals of motion is the quantum transfer matrix T^'(u), which depends on the spectral 
parameter m£C and the Young diagram A. It is obtained as the (super)trace of the quantum monodromy 
matrix (u) in the auxiliary space carrying the irreducible representation of the symmetry algebra labeled 
by A: 

T< A >(u)=str aux T< A >( W ). 

We deal with covariant tensor irreducible representations (irrcps) of the supcralgebra. 

The transfer matrices for different A are known to be functionally dependent. For rectangular diagrams 
A = (s a ) with a rows and s columns the functional relation takes the form of the famous Hirota difference 
equation: 

T(a, s, u+l)T(a, s, u— 1) — T(a, s+1, u)T(a, s—l,u) = T(a+l,s, u)T{a — 1, s, u), (1-1) 

where the transfer matrix for rectangular diagrams, after some a, s-dependent shift of the spectral parameter 
u, is denoted by T(a,s,u). We see that it enters the equation as the r-function [29] . Since all the T's commute, 
the same relation holds for their eigenvalues. We use the normalization in which all non-vanishing T's are 
polynomials in u of degree N, where N is the number of sites in the spin chain. 

While the functional relation is the same for ordinary and super-algebras, the boundary conditions are 
different. In the gl(K)-c&se, the domain of non-vanishing T's in the (s,a) plane is the half-strip s > 0, 
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< a < K. In the gl(K\M)-ca.se, the domain of non-vanishing T's has the form of a "fat hook". It is shown 
in Fig. 1. To ensure compatibility with the Hirota equation, the boundary values at s = and a = should 
be rather special. In our normalization, 

T(0, s, u) = 4>{u - s) , T(a, 0, u) = (f>(u + a) , (1.2) 

where 4>{u) = Y[i( u ~ ®i) IS a fixed polynomial of degree N which characterizes the spin chain. 

Our goal is to solve the Hirota equation, with the boundary conditions given above, using the classical 
methods of the theory of solitons. This program for the case of ordinary Lie algebras (and for models with 
elliptic i?-matrices) was realized in [16]. In this paper, we extend it to the case of supcralgcbras (for models 
with rational i?-matrices). 

One of the key features of soliton equations is the existence of (auto) Backlund transformations (BT), 
i.e., transformations that send any solution of a soliton equation to another solution of the same equation. A 
systematic way to construct such transformations is provided by considering an ovcr-dctcrmincd system of 
linear problems (called auxiliary linear problems) whose compatibility condition is the non-linear equation at 
hand. We introduce two Backlund transformations for the Hirota equation. They send any solution with the 
boundary conditions described above to a solution of the same Hirota equation with boundary conditions of 
the same fat hook type but with different K or M. Specifically, one of them lowers if by 1 and the other one 
lowers M by 1 leaving all other boundaries intact. We call these transformations BTl and BT2. Applying 
them successively K + M times, one comes to a collapsed domain which is a union of two lines, the s-axis 
and the a-semi-axis shown in Fig.l, meaning that the original problem gets "undressed" to a trivial one. 
This procedure appears to be equivalent to the nested Bethe ansatz. Different orders in which we diminish 
K and M give raise to different "dual" systems of nested Bethe Ansatz equations. All of them describe the 
same system. They correspond to different choices of the basis of simple roots. 

Let k, m be indices running from to if and from to M respectively, and let Tfc jrn (a, s, u) be the transfer 
matrices obeying the Hirota equation with the boundary conditions as above but with if, M replaced by 
k, m. They are obtained from Tk,m{<Ii s, u) by a chain of BT's. Namely, Tfc_i jm 's and T^m—i's are solutions 
to the auxiliary linear problems for the Hirota equation for Tfc r „'s: 

rp BT2 „ 



BTl 



T, 



k— l,m 

The explicit formulas of these transformations are given below (eq.(3.29) and eq.(3.30)). The functions 
Tk, m {o-: s, u) are polynomials in u for any a, s, k, m but the degree depends only on k, m. Let Qk,m(u) be the 
boundary values of the T^ m {a, s, u), i.e., 

T fcim (0, s, u) = Qk, m (u - s) , T kiTn (a, 0, it) = Q k , m (u + a) . 

However, they are not fixed for the values of k,m other than k = K , m = M (when Qk,m( u ) = an d 
k = m = (when Qo i o(w) is put equal to 1) but are to be determined from a solution to the hierarchy of 
Hirota equations. In fact these Q's are Baxter polynomial functions whose roots obey the Bethe equations. 

The result of a successive application of the transformations BTl and BT2 does not depend on their 
order. This fact can be reformulated as a discrete zero curvature condition 



U k,m+i V k,m - Vk+i,mU k rn (1.3) 



for the shift operators in k and m: 



jj^ / \ _ Qk+l,m(u) Qk,m\U+2) ^ 2 a„ 
' m Qk+l,m(u + 2)Q k , m (u) 

?/ QkM u ) Qfc,m+i(M+2) 2du 

V k , m {u) = -= -. s n 7-T - e . (1.4) 
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Figure 2: The R-matrix. 



Relation (1.3) is equivalent to the following Hirota equation for the Baxter Q-functions: 

Qk,m(u)Qk+l,m+l{u + 2) - Qk+l,m+l(u)Qk,m(u + 2) = Q k,m+l(u)Q k+l,m(u + 2). (1.5) 

This equation represents our principal new result. By analogy with Baxter's TQ-relations, we call eq. (1.5) 
the QQ-relation. It provides the most transparent way to derive different systems of Bcthc equations for the 
generalized spin chain and "duality transformations" between them. We also show that a number of more 
general Hirota equations of the similar type (i.e., acting in the space spanned by k, m and a particular linear 
combination of a,s,u) hold for the full set of functions 2\ m (a, s, u). They lead to a system of algebraic 
equations for their roots which generalizes the system of Bcthc equations. 

The transfer matrices can be expressed through the Q-functions via generalized Baxter's TQ-relations. A 
simple way to represent them is to consider the (non-commutative) generating scries of the transfer matrices 
for one-row diagrams, 

oo 

W k , m {u) =Q^ m {u)Y,T k , m (l,s,u+s + l)e 2sd ", (1.6) 

where the factor in front of the sum is put here for the proper normalization. The operator W is similar 
to the wave (or dressing) operator in the Toda lattice theory. We prove the following fundamental operator 
relations: 

W k+ltm = UkXWk,m , W k , m+1 = V Km W Km , (1.7) 

which implement the Backlund transformations on the level of the generating series. These relations allow 
one to represent Wk,m (the quantity of prime interest) as an ordered product of the operators Uu m and 

Vk,m along a zigzag path from the point (0, 0) to the point (K, M) on the (k, m) lattice. This representation 
provides a concise form of the generalized Baxter relations. Different zigzag paths correspond to different 
choices of the basis of simple roots, i.e., to different "dual" forms of supersymmetric Bethe equations (see 
Fig. 17). 

The normalization (1.2), with the roots 6i of the polynomial 4>{u) being in general position, implies that 
the "spins" of the inhomogeneous spin chain belong to the vector representation of the algebra gl(K\M). 
Higher representations can be obtained by fusing of several copies of the vector ones. According to the 
fusion procedure, the corresponding #'s must be chosen in a specific "string-like" way which means that 
the differences between them are even integers constrained also by some more specific requirements. This 
amounts to the fact that the Q-functions become divisible by certain polynomials with explicitly prescribed 
roots located according to a similar string-like pattern. If one redefines the Q-functions dividing them by 
these normalization factors, then the QQ-relation gets modified and leads, in the same way as before, to the 
systems of Bcthe equations with non-trivial right hand sides (sometimes called "vacuum parts"). 

In section 7, we specify our construction for two popular examples, the gZ(l|l) and <?Z(2|1) superalgebras. 
We present the TQ- and QQ-rclations for all possible types of Q-functions Q km (u) (0 < k < K, < m < 
M) as well as Bethe equations for all types of simple roots systems and representations of the superalgebra 
(including typical representations with a continuous component of the Kac-Dynkin label) and compare them 
with the results existing in the literature. In our formalism, the construction becomes rather transparent and 
algorithmic. In this respect, our method is an interesting alternative to the algebraic Bethe ansatz approach. 
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Figure 3: The Yang-Baxter equation. 

2 Fusion relations for transfer matrices and Hirota equation 



2.1 Quantum transfer matrices 

Let V = C K © C M = C K \ M be the graded linear space of the vector representation of the superalgebra 
gl(K\M). The fundamental gZ(A"|M)-invariant rational i?-matrix acts in V <g> V and has the form [1, 2] 



R(u) = ul + 277 = ul + 2 



K+M 

E 

7,7'=1 



(-1) 



P(7 ) 



77 



7 7 • 



(2.1) 



Here 7 is the identity operator and 77 is the super-permutation, i.e., the operator such that U(x ® y) = 
(—iy(x)p(y)y (g, x f or an y homogeneous vectors x,y £ V, and E a p are the generators of the (super) algebra. 
In components, they read 



HP' 



\Ey*yt ) Q / — 3a^0~a' 



Note that 77 2 = I, as in the ordinary case. The notation p(x) is used to denote parity of the object x 
(see Appendix A). The R(u) has only even matrix elements. The complex variable u is called the spectral 
parameter. The R- matrix obeys the graded Yang-Baxter equation (Fig. 3): 



J2 (-iy>MW'MT' , »R#(u - v)R}«';,(u)R},};,(v) 

77'7" 



E 

77'7" 



(_i)P(70(p(/3^+ P (7'0) i? ^; (w )^7;; (M) ^77; (u _ v) 



(2.2) 



This is the key relation to construct a family of commuting operators, as is outlined below. 

The i?-matrix is an operator in the tensor product of two linear spaces (not necessarily isomorphic). It 
is customary to call one of them quantum and the other one auxiliary space. In Fig. 2 they are associated 
with the vertical and the horizontal lines respectively. Let us fix a set of N complex numbers 6i and consider 
a chain of N fundamental i?-matrices R(u — &i) with the common auxiliary space V. Multiplying them 
as linear operators in V along the chain, we get an operator in V® N which is called quantum monodromy 
matrix. In components, it can be written as 



T 70 '<"\ } ( u ) 

7w,{/3i} v ' 



E 



p7N-iajv , 



(u-9 N ) . ..R^u-e^R^iu-ex) (_i)5^WaO+p(ft)]I&(ai) (2.3) 



It is usually regarded as an operator-valued matrix in the auxiliary space (with indices 70, 7at), with the 
matrix elements being operators in the quantum space V® N (see Fig. 4). Setting 70 = jn = 7 and summing 
with the sign factor (— 1)p(t), we g e ^ ^ ne supertrace of the quantum monodromy matrix in the auxiliary 
space, 

T(u) = strT(u) = ^(-l) p(7) 77(u) (2.4) 

7 



7 



c/ YN ,{Pi} - Y ° 
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^ Yi 


'p. 









Yn-1 



Yn 



01 Q 2 c N 

Figure 4: TVie monodromy matrix. 



which is an operator in the quantum space (indices of the quantum space arc omitted). It is called the 
quantum transfer matrix. Using the graded Yang-Baxter equation, it can be shown that the T(u) commute 
at different u: \T(u), T(u')] = 0. Their diagonalization is the subject of one or another version of Bethe 
ansatz. 

The (inhomogeneous) integrable "spin chain" or a vertex model on the square lattice is characterized by 
the symmetry algebra gl(K\M) and the parameters 9i (inhomogeneities at the sites or "rapidities"). Given 
such a model, the family of transfer matrices T{u) (2.4) can be included into a larger family of commuting 
operators. It is constructed by means of the fusion procedure. 



2.2 Fusion procedure 

Using the 7?-matrix (2.1) as a building block, it is possible to generate other gl(K\M )-invariant solutions to 
the (graded) Yang-Baxter equation. They are operators in tensor products of two arbitrary finite dimensional 
irreps of the symmetry algebra. We consider covariant tensor irreps which are labeled (although in a non- 
unique way, see Appendix A) by Young diagrams A. Let V\ be the corresponding representation space, then 
the 7?-matrix acts in V\ ® Vy ■ The construction of this 7?-matrix from the elementary one is referred to as 
fusion procedure [7, 32], [33]- [36]. 

Here we consider fusion in the auxiliary space, which allows one to construct, by taking (super)trace 
in this space, a large family of transfer matrices commuting with T(u). In words, the fusion procedure 
consists in tensor multiplying several copies of the fundamental 7?-matrix (2.1) in the auxiliary space and 
subsequent projection onto an irrep of the gl[K\M) algebra. As a result, one obtains an 7?-matrix RS x \u) 
whose quantum space is still V (the space of the vector representation of gl(K\M)) and whose auxiliary 
space is the space V\ of a higher irrep A of this algebra. The fact that the 7?-matrix obtained in this way 
obeys the (graded) Yang-Baxter equation follows from the observation [7, 34] that the projectors onto higher 
irreps can be realized as products of the fundamental 7?-matrices (2.1) taken at special values of the spectral 
parameter. These values correspond to the degeneration points of the 7?-matrices. 

In particular, the 7?-matrix (2.1) has two degeneration points u = ±2. Indeed, one can easily see that 
det R(u) = {u + 2) d +(u-2) d - , d± = ^((K + Mf ± (K - M)) , (2.5) 

R{±2) = ±2(7 ± 77) = ±4P± , (2.6) 

where P± are projectors onto the symmetric and antisymmetric subspaces in V ®V . The dimensions of 
these subspaces are equal to d±, and the projectors are complimentary, i.e., P+P- — P-P+ — 0. Therefore, 
to get the 7?-matrix with the auxiliary space carrying the A = (2 1 ) irrep (respectively, the A = (l 2 ) irrep), 
one should take the tensor product R(u + 2) ® R(u) (respectively, R(u — 2) ® R(u)) in the auxiliary space 
and apply the projector P + (respectively, P-): 

R (2l) (u) =P+\r(u + 2)®R(u)}p + , R {l2) (u) = P_ \r(u - 2) ® R(u)] P- . 
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Figure 5: The Young diagram decorated by spectral parameters. 



Here, the tensor product notation still implies the usual matrix product in the quantum space. Note that 
fj(2 )(— 2) vanishes identically since the projector P+ gets multiplied by the complementary projector P_ 
coming from one of the i?-matrices in the tensor product. Similarly, R^ 1 '(2) = 0. 

In a more general case, the procedure is as follows. Let A be the Young diagram associated with a given 
irrep of the algebra 1 . Let n — |A| be the number of boxes in the diagram and let 

be the projector onto the space of the irreducible representation. Write in the box with coordinates (i-th 
line counting from top to bottom and j'-th column counting from left to right) the number — u~ 2(i — j), 
as is shown in Fig. 5. Then R( x \u) is constructed as 

R^(u) =P A [ (g) i2(uy)] Px (2.7) 
(y)eA 

Here, the tensor factors are placed from right to left in the lexicographical order, i.e., elements of the first 
row from left to right, then elements of the second row from left to right, etc. For example, for the diagram 
A = (3, 2) the product under the projectors reads R(u) ® R(u — 2) ® R(u + 4) ® R(u + 2) <E> R(u). 

Omitting further details of the fusion procedure (see, e.g., [35, 36]), we now describe the structure of zeros 
of the fused i?-matrix R( x \u) which is essential for what follows. As we have seen, the i?-matrix R( x \u) 
is obtained by fusing n — |A| fundamental R- matrices in the auxiliary space. Then, by the construction, 
any matrix element of the R^(u) is a polynomial in u of degree < n. However, one can see that at some 
special values of the spectral parameter the complimentary projectors get multiplied and the matrix R^(u) 
vanishes as a whole thing. In other words, the matrix elements have a number of common zeros. We call 
them "trivial zeros" . A more detailed analysis shows that the number of trivial zeros is n — 1 (counting with 
multiplicities) and thus a scalar polynomial of degree n — 1 can be factored out. Namely, this polynomial is 
equal to the product of the factors Uij — u — 2(i — j) over the boxes of the diagram A with (i, j) G A except 
i=j = 1. 

We consider here only the representations of the type A = (s a ) corresponding to rectangular Young 
diagrams with s rows and a columns (rectangular irreps). In this case 

R {sa} (u)^r a , s (u)R^\u), 

where matrix elements of R^ s '(u) are polynomials of degree < 1 and 

s a 

r a , s (u) = u- 1 Y[Y[(u + 2j -21), a,s>l, (2.8) 

3 = 11=1 



1 For superalgcbras this may be a delicate point since this correspondence is in general not one-to-one. 
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is the polynomial of degree as — 1. For a future reference, we give its representation through the Barnes 
function G{z) (a unique entire function such that G(z + 1) = T(z)G(z) and G(l) = 1): 

r»W = 2-„- g'?^'^ G f-°^ . (2.9) 
G(±u + 1) G(|u + s — a + 1) 

The quantum monodromy matrix with the auxiliary space V\ is defined by the same formula as 
before but with the i?-matrix RS X '. The transfer matrix is 

T^\u) = stT Vx T^{u) (2.10) 

The transfer matrices commute for any A and u: [T^ x ' (u), > (it')] = 0. This commuting family of operators 
extends the family (2.4) which corresponds to the one-box diagram A. For the empty diagram A = we 
formally put TW(tt) = 1. For rectangular diagrams A = s a we use the special notation T^ x \u) := T*(u). 



2.3 Functional relations for transfer matrices 

The transfer matrices are functionally dependent. It appears that all T^\u) can be expressed through 
Tg(u) or Ti(u) by means of the nice determinant formulas due to Bazhanov and Reshetikhin [9]. They are 
the same for all (super) algebras of the type gl(K\M) with K,M > [12]: 



T^fu) = det Tl , , Au-2i + 2) 



det T^ +l 3 (u + 2i-2). 

l<i,i<Ai 



(2.11) 



For rectangular diagrams, they read: 



T«(u) = det Tl H _ j (u-2i + 2) 



(2.12) 

= det T? +i - J (u + 2i-2). 

These formulas should be supplemented by the "boundary conditions" T s °(u) = Tq(u) = 1, T~ n (u) = 
T°: n (u) = for all a, s > and n > 1. Since all the transfer matrices commute, they can be diagonalized 
simultaneously, and the same relations are valid for any of their eigenvalues. Keeping this in mind, we will 
often refer to the transfer matrices as scalar functions and call them "T-functions" . 

Applying the Jacobi identity for determinants to formulas (2.12) (see, e.g., [37, 16]), one obtains a 
closed functional relation between transfer matrices for rectangular diagrams, equivalent to Hirota difference 
equation: 

T«(u - 2)T S » - T« +1 (u - 2)2?_ 1 (u) = T^ x (u - 2)T«+\u) (2.13) 

The bilinear form of the functional relations was discussed in [10], [11], [16], [37]- [40]. Below we illustrate 
this relation by simple examples. 



2.3.1 Simple examples of the functional equation 

To figure out the general pattern, it is useful to consider the simplest case of the fusion of two transfer 
matrices with fundamental representations in the auxiliary space. The proof is summarized in Fig. 6. Let 
us represent the first term in eq.(2.13) as 

Tl(u - 2)Tl{u) = tr aux , vw Pao(u - 2)7i (u)] (2.14) 

The index denotes the common quantum space represented in the figure by several vertical lines, the indices 
1, 2 denote the two copies of the auxiliary space V. Inserting I = P + +P_ inside the trace, we get two terms. 
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Figure 6: Fusion of two fundamental transfer matrices and the corresponding Hirota equation. As usual, each 
line carries the vector representation, each crossing corresponds to the insertion of the R-matrix, and the 
closed lines correspond to taking traces in the auxiliary space. The insertions of projectors P^ are interpreted 
as insertions of the R-matrices with given values of spectral parameters. 



Using the projector property P± = P±, we immediately see that the term with P + is equal to T 2 1 (u — 2), by 
the definition of the latter. The term with P_ does not literally coincide with the definition of T^(u) since 
the order of the horizontal lines is reversed. However, plugging _P_ = — jR(— 2) and using the Yang-Baxter 
equation to move the vertical lines to the other side of the i?-matrix 2). we come to the equivalent graph 
with the required order of the horizontal lines. Finally, we get 

Tl(u - 2)T 1 1 (m) = T\(u - 2) + Tl(u) . (2.15) 

Hence we have reproduced the simplest case of eq.(2.13). 

Let us comment on the general case. An illustrative example is given in fig. 7. The rectangular Young 
diagrams are decorated by the values of the spectral parameter 

Uij =u-2(i-j), (2.16) 

as is explained above. Each box of such decorated diagrams corresponds to a line characterized by a given 
value of the spectral parameter. All these lines cross each other and the i?-matrices are associated to the 
crossing points. The order of the crossings is irrelevant due to the Yang-Baxter equation. The reshuffling of 
spectral parameters in the second and third terms of equation (2.13) (corresponding to the second and the 
third lines in the figure) is done by the exchange of a row or a column between the two diagrams of the first 
line. The decoration of the resulting Young diagrams follows the rule (2.16). Of course this is not a proof, 
just an illustration. 



2.3.2 Analogy with formulas for characters and identities for symmetric functions 

Equations (2.11) are spectral parameter dependent versions of the second Weyl formula for characters of 
(super)groups (see [20]): 

X\(wi, ■ ■ ■ , wk;vi, . . . ,vm) = det S\ i+ i-j . ( 2 -17) 

i,j=l,...,K+M 

In the theory of symmetric functions such formulas are known as the Jacobi-Trudi determinant identities. 
Here Aj are the lengths of the rows of the diagram A, S n are super Schur polynomials defined as 



n m =l (1 - ZV m) = 



u k n , (2-18) 

Uk=l l 1 - zw k) 



and Wk, i) m are eigenvalues of the A- and Z?-parts of the diagonalized element of the supergroup in the 
matrix realization of the type (Al). 
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Figure 7: The Hirota-type relation illustrated by rectangular Young diagrams decorated by spectral parameters. 
To obtain the first term in the r.h.s., one takes the first row of the second diagram in the l.h.s. and puts it 
on the top of the first one. To obtain the second term in the r.h.s., one takes the first column of the first 
diagram in the l.h.s. and attaches it to the second one from the left. 



det Ss+i-j- 

i,j=l,...,a 



For the rectangular irrep with A = s a eq. (2.17) reads 

Xs° = x(a,s) 

The characters Xs a satisfy the bilinear relation 

X 2 (a, s ) = X( a + 1) s )x( a ~ 1, s) + xi a , s + 1)x(Oj s - !) 



(2.19) 



(2.20) 



which follows from the Jacobi identity for determinants. It is the spectral parameter independent version of 
the Hirota equation. 



2.3.3 Normalization and boundary conditions 

A simple redefinition of the T-functions by a shift of the spectral parameter, 

T^ s \u)=T^(u-s + a) (2.21) 

brings equation (2.13) to the form 

T^ s \u + l)T^ s \u - 1) - T( a ' a+1 \u)T<- a >'- l \u) = T^+^iu^-^iu) (2.22) 

which appears to be completely symmetric with respect to interchanging of a and s. This is the famous 
Hirota bilinear difference equation [18] which is the starting point of our approach to the Bethe ansatz and 
generalized Baxter equations in this paper. For convenience, we also give the determinant representation 
(2.12) in terms of T< Q ' s )(w): 



T M( U ) = det T^' s+i - j \u + a+l-i-j) 

l<i,j<a 

det T {a+l - 3A \u + s + l-i-j). 

X<A,j<s 



(2.23) 



Let us comment on the meaning of equation (2.22). On the first glance, there is no much content in this 
equation. Its general solution (with the boundary conditions fixed above) is just given by formulas (2.23) 
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with arbitrary functions T^ 1,s ^(u) or T( a,1 )(it). However, in the problem of interest these functions are by no 
means arbitrary. They are to be found from certain analytic conditions. For the finite spin chains with finite 
dimensional representations at each site these conditions simply mean that T^- a,s \u) must be a polynomial 
of degree asN, where N is the length of the chain, with (as — 1)N fixed zeros. These zeros are just the trivial 
zeros coming from the fusion procedure. Their location is determined by the scalar factor r ajS (it) (2.8) of 
the i?-matrix. We thus see that the polynomial T^ a,s '(u) for all a, s > 1 must be divisible by the polynomial 

N 

~[r a . s (u -8 + a-di) := $(a,s,u) . 

i=l 

This constraint makes the problem non-trivial. 
Let us introduce the function 

N 

0(«) = IK" - ^ ' ( 2 ' 24 ) 

i=l 

in terms of which the scalar polynomial factor is written as 

N s a 

$(a, s, u) = r a , s (u - s + a - 6i) = (f>~ l (u - s + a) T ~~\<f>(u - s + a + 2j — 21) , a, s > 1. 
i=l 3=1 1=1 

The representation through the Barnes function (2.9) allows us to extend this formula to all values of a and 

K ' V Li G(%(u-s + a) + l-6i) G(%(u + a-a) + l-6i) K ' 

Extracting it from the T^(u), we introduce the T-function 

T(a,s,u) = <f>-\a,s,u)T ia - s) (u) (2.26) 

which is a polynomial in u of degree N for all a, s > 0. Note that at a — or s — equation (2.25) yields 
$(0, s, u) = l/<p(u - s), $(0,0, u) = l/4>{u + a), so 

T(0,s,u)=<f>(u-s), T(a,0,u) = (j)(u + a). (2.27) 

It is important to note that the rcnormalizcd T-function T(a, s, u) obeys the same Hirota equation (2.22) 
as the T-function T^ a ' s ^(u). Indeed, it easy to check that the transformation 

_^ foiu + s + ^hiu + s-^hiu-s + ^Uu-s-^T^iu), (2.28) 

where fi are arbitrary functions, leaves the form of the equation unchanged. Equation ((2.25)) shows that 
the function $(a, s, u) is precisely of this form (the factor 2 asN is easily seen to be of this form, too). 

The main difference between the "bosonic" gl(K) and supersymmetric gl(K\M) cases is in the boundary 
conditions for the transfer matrices in the (s, a)-plane. For the algebra gl(K) the rectangular Young diagrams 
live in the half-band s>0,Q<a<K, while for the superalgebra the rectangular diagrams live in the domain 
shown below in Fig. 9. The Hirota equation with boundary conditions of this type will be our starting point 
for the analysis of the inhomogeneous quantum integrable super spin chains and it will allow us to obtain 
the full hierarchy of Baxter relations, the new Hirota equation for the Baxter functions and the nested Bcthe 
ansatz equations for all possible choices of simple root systems for gl(K\M). This naturally generalizes the 
known relations for the spin chains based on the gl{K) algebra. 



3 Hierarchy of Hirota equations 

Throughout the rest of the paper we deal with rectangular irreps only and use the normalization (2.26), 
where all the "trivial" zeros of the transfer matrix are excluded. This normalization was used in [16, 17], for 
the supersymmetric case sec [15]. 
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Figure 8: Hirota equation in the (s, a)-plane. 



3.1 Hirota equation and boundary conditions for superalgebra gl(K\M) 

As we have seen in the previous section, the functional relations for transfer matrices of integrable quantum 
spin chains with "spins" belonging to representations of the superalgebra gl(K\M) can be written [12] in 
the form of the same Hirota equation as in the case of the ordinary Lie algebra gl{K): 

T(a, s, u+l)T(a, s, u— 1) — T(a, s+1, u)T(a, s—l,u) = T(a+1, s, u)T(a — \, s,u). (3-1) 

It can be schematically drawn in the (s, a)-plane as is shown in Fig. 8. All the non- vanishing T(a, s, u)'s are 
polynomials in u of one and the same degree N equal to the number of sites in the spin chain. 

To distinguish solutions relevant to Bethe ansatz for a given (super)algcbra, one should specify the 
boundary conditions in the discrete variables a and s. For the superalgebra gl(K\M) one can see [41, 12] 
that 

T(a,s,u) = if : (i) a < or (ii) s < and a + 1 0, or (iii) a > K and s > M. (3.2) 

(see Fig. 9, where we use the letters k, m rather than K, M for later references). The latter requirement, that 
T(a,s,u) vanishes if simultaneously a> K and s > M, comes from the fact that the Young superdiagrams 
for gl(K\M) containing a rectangular subdiagram with K + 1 rows and M + 1 columns are illegal, i.e., the 
corresponding representations vanish [21]. Note that we want the Hirota equation to be valid in the whole 
(s, a) plane, not just in the quadrant a, s > 1. This is why we have to require that T(0, s, u) docs not vanish 
identically on the negative s-axis, otherwise the Hirota equation would break down at the origin a = s = 0. 

The boundary values of T(a, s, u) are rather special. For example, at a = eq. (3.1) converts into 

T(0, s, ti+l)T(0, s, u-l) = T(0, s + 1, u)T(0,s-l, u) 

which is a discrete version of the d'Alembert equation with the general solution T(0, s, u) = f+(u+s)f—(u—s) 
where f± are arbitrary functions. In the normalization (2.26) we have f+(u) = 1, f-(u) = <fr{u) (see (2.27)). 
Similarly, the boundary function at the half-axis s = 0, a > is normalized to depend on u + a only in which 
case it has to be equal to f-{u + a) = (p(u + a). As soon as this is fixed, there is no more freedom left, 
and the boundary functions at the interior boundaries 2 are in general products of a function of u + s and a 
function of u — s on the horizontal line (respectively, of u + a and u-aon the vertical line) . One more thing 
to be taken into account is the identification (up to a sign) of the T-functions on the two interior boundaries: 

T(K, M + n,u) = {—l) nM T(K + n, M, u), n > 0. (3.3) 

This equality reflects the fact that the two rectangular Young diagrams of the shapes ((M+n) K ) and 
(M K+n ) correspond to the same representation of the algebra gl(K\M) with the Kac-Dynkin label b\ = 
. . . = bic-i = 0, &a" = M + n, bx+i = . . . = bic+M-i = [21]. Note also that every point with integer 
coordinates inside the domain in Fig. 9 corresponds to an atypical representation while the points on the 
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Figure 9: Boundary conditions for the Hirota equation in the case of the superalgebra gl(k\m). Note that 
T(0,s,u) = Qk, m (u — s) 07i the whole axis — oo < s < +oo. 



interior boundaries correspond to typical ones, and in this respect the coordinate along the boundary can 
be treated as a continuous number. 



Summarizing, we can write 

T(0,s,u) = 
T(a,0,u) 
T(K,s,u) 
T{a, M, u) 



Qk,m(u — s), — oo < s < oo, 
Qk,m(u + a), < a < oo, 

Qk,q(u + s + K)Q ,m(u - s - K), M < s < oo, 
{-l) M{a - K) Q K ,o{u + a + M)Q ,m(u -a- M), K<a< oo, 



(3.4) 



where the polynomial boundary function Qk,m(u) — 4>(u) is regarded as a fixed input characterizing the 
quantum space of the spin chain while the polynomials Qk,o(u) and Qo,m( u ) are to be determined from the 
solution to the Hirota equation. At K — M — the domain of non-vanishing T's shrinks to the axis a = 
and the half-axis s = 0, a > 0, and the "gauge" freedom allows one to put Qo,o( w ) = 1- It should be noted 
that the identification (3.3) does not yet imply the coincidence of the Q-functions in the third and the fourth 
lines of eq. (3.4) (i.e., on the horizontal and vertical parts of the interior boundary). In fact the specific 
form of the boundary conditions given in eq. (3.4) (as well as the sign factor (— l)^* 1- -^)) is determined 
by a consistency with a more general hierarchy of Hirota equations connecting the T-functions for different 
values of K and M. The uniqueness of the boundary conditions (3.4) will be justified in the next subsection. 

In the case of the usual Lie algebra gl{K) (M = 0) the boundary conditions (3.4) become the same as the 
ones imposed in [17] (in the original paper [16] a gauge equivalent version was used). The domain of non- 
vanishing T's in Fig. 9 degenerates so that the vertical strip collapses to a line. Therefore, the T-functions 
on the interior boundaries cease to be dynamical variables and become equal, up to a shift of the argument, 
to the fixed function Qk.q = </>(w) characterizing the spin chain (see [16, 17]). 

One of the possible setups of the problem is the following: we fix the polynomial Qk,m( u ) an d then solve 
the Hirota equation with the aforementioned boundary and analytic conditions. The result is a finite set of 
solutions for T(a, s,u) which yield the spectrum of eigenvalues of the quantum transfer matrix. We proceed 
by constructing a hierarchy of Hirota equations connecting neighboring "levels" of the array in Fig. 9, i.e., 
equations connecting T-functions for which K or M differ by 1. The existence of such a hierarchy follows 
from classical integrability of the Hirota equation. Decreasing K and M by 1, one can "undress" step by 
step the original gl(K\M) problem to an empty problem formally corresponding to g/(0|0). This procedure 
appears to be equivalent to the hierarchial (nested) Bcthe ansatz. 



2 We call the boundaries at a 
interior ones. 



and s = 0, a > exterior and the boundaries inside the right upper quadrant in Fig. 9 
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3.2 Auxiliary linear problems and Backlund transformations 



Like almost all known nonlinear intcgrable equations, the Hirota equation (3.1) serves as a compatibility 
condition for over-determined linear problems [30, 16]. To introduce them, it is convenient to pass to the 
new variables 

p = |(u — s — a) 



q = |(u + s + a) 



(3.5) 



r = h(—u — s + a) 

We call them "chiral" or "light-cone" variables while the original ones will be refered to as "laboratory" 
variables. Here are the formulas for the inverse transformation, 

a = q + r , s = —p — r , u = p + q (3-6) 

and for the transformation of the vector fields: 

d p = d u -d s , d q = d u + d a , d r = d a - d s (3.7) 

We set r(p, q, r) = T(q + r, — p — r, p + q) and introduce the following linear problems for an auxiliary 
function ip = ip(p, q, r): 

d r , r(p+l,r + l) r 



r(p+l)r(r + l) 
r(g+l,r+l) r 



V> = ip(p + i) 



(3.8) 



r(q + l)r(r + 1), 

where we indicate explicitly only those variables that are subject to shifts. The compatibility means that 
the difference operators 



r(p+l, r+1) r 



and 



,9 r _ r(g + l,r+l) r 
r(g + l)r(r + l) 



r(p + l)r(r+ 1), 

commute. This leads to the relation 

r(p + l)r(q + 1, r + 1) + r(q + l)r(p + 1, r + 1) = ft(2p, 2g)r(r + l)r(p +l,q+l) 

where h can be an arbitrary function of p and q. In the original variables this equation reads 

T(a + l)T(a - 1) + T(s + l)T(s - 1) = h(u-s-a, u+s+a) T(u+ l)T(u - 1) 

From the boundary conditions (3.4) at a = or s = it follows that h = 1 and we obtain the Hirota equation 
(3.1). 

An advantage of the "light-cone" variables is their separation in the linear problems: the first problem 
does not involve q while the second one does not involve p. However, in contrast to the "laboratory" 
variables a, s,u, they have no immediate physical meaning. Coming back to the "laboratory" variables, we 
set ip(Py 9; r ) = ^(<Z + r i ~ P ~ r ; V + o) an d rewrite the linear problems (3.8) in the form 



oda _ dB r(a-l,s + l,M)r(q,s-l,M+l) 
T(a, s, u) T(a— 1, s, u+1) 

= g a _a s T(o— 1, s+1, u)T(a+l, s, u+1) 
T(a, s, u) T(a, s + 1, u + 1) 



*(a-l,s + l,u) = *(a— l,s,u+l) 



*(a-l, s + 1, it) = *(a, s + 1, u + 1) 



(3.9) 



Because the T-functions can vanish identically at some a, s, we eliminate the denominators by passing to 
the new auxiliary function F = Tty, in terms of which we have 



T{a — 1, s, u + l)F(a, s,u) + T(a, s — 1, u + l).F(a — 1, s + 1, u) = T(a, s, 14)^(0— 1, s, u+1) 
T(a, s + 1, w + l)F(a, s, u) — T(a + 1, s, 14+1)^(0— 1, s+1, u) = T(a, s, u)F(a, s + 1, u+1) 



(3.10) 
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Figure 10: Boundary conditions for the first Bdcklund transformation in the (s, a) plane. 



Note that the second equation can be obtained from the first one by the transformation T(a, s, u) — ► 
( — l) as T(— s, —a, u) (and the same for F) which leaves the Hirota equation invariant. Note also that the pair 
of equations can be written in a matrix form as follows: 



T(a — l,s,u) T(a,s — l,u) 
T(a, s + l,u) —T(a+l,s,u 




F(a, s, u — 1) 
F(a-l,s+l,u-l) 



T(a, s,u—l) 



F(a—1, s, u) 
F(a, s+1, u) 



(3.11) 



A remark on the symmetry properties of the linear problems is in order. One can see that while the 
Hirota equation (3.1) written for the function (— l)2( a +s ~)T(a,s,u) is form-invariant with respect to any 
permutation of the variables a, s, u and changing sign of any variable, the system of two linear problems 
(3.10) is not. To make the symmetry explicit, we multiply both sides of eq. (3.11) by the matrix inverse to 
the one in the left hand side and use the Hirota equation for T's. In this way we get another pair of linear 
problems, 



T(a+1, s+l, u)F(a, s, it) — T(a+1, s, u + l)F(a, s+1, u— 1) = T(a, s, u)F{a+l, s+1, u) 
T(a, s,u+l)F(a, s,u — l) — T(a, s — 1, u) F(a, s + 1, u) = T(a+1, s,u) F(a—1, s,u) 



(3.12) 



which are equivalent to (and thus compatible with) the pair (3.10) by construction. The set of four linear 
problems (3.11), (3.12) possesses the required symmetry (note that the second equation in (3.12) is symmetric 
by itself, and its structure resembles the Hirota equation). Furthermore, the Hirota equation can be derived 
as a compatibility condition for any two linear problems of these four, and the other two hold automatically. 
The four linear problems (3.11), (3.12) can be combined into a single matrix equation: 



/ o 

—T(a, s, it — 1) 
T(a, s+1, u) 
\ —T(a+l,s,u) 



T(a,s,u—1) —T(a,s+l,u) T(a + l,s,u) \ 

T(a—l,s,u) T(a,s — l,u) 

—T(a—l,s,u) — T(a,s,u+1) 

—T(a,s—l,u) T(a,s,u+1) J 



( F(a-l,s,u) \ 
F{a, s + 1, u) 
F{a, s, u—1) 
\ F(a-l,s + l,u-l) J 



= 0. 



(3.13) 

The Hirota equation implies that the determinant of the antisymmetric matrix in the left hand side vanishes. 
If this holds, the rank of this matrix equals 2, so there are two linearly independent solutions to the linear 
problem (3.13). Their meaning will be clarified below. The symmetric form of the linear problems was 
suggested in [42] . For more details on the linear problems for the Hirota equation and their symmetries see 
[30, 31, 42, 43] and Appendix B. 

There is a remarkable duality between T(a,s,u) and F(a,s,u) [30, 16]: one can exchange the roles of 
the functions T, F and treat eqs. (3.10) as an over-determined system of linear problems for the function T 
with coefficients F. Their compatibility condition is the same Hirota equation for F: 



F(a, s,u+ l)F(a, s,u — 1) — F(a, s + 1, u)F(a, s — 1, u) = F(a + 1, s, u)F(a — 1, s, u) 



(3.14) 
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Figure 11: Two equations for the first Backlund transformation (BT1) in the (s,a) plane: the variables 
change within a "brick" consisting of two neighboring squares. 



This can be seen by rewriting (3.10) in yet another equivalent form. Namely, shifting a — > a + 1 in the first 
equation and s — > s — 1 in the second, we represent the two equations in the matrix form 



F(a+l,s,u) F(a,8+l,u) 
F(a,s — l,u) —F(a—l,s,u) 



T(a, s, it+1) 
T(a+l,s-l,u+l) 



F(a, s, u+1) 



T(a+1, s, u) 
F(a, s — l,u) 



(3.15) 



which is to be compared with (3.11). It is obvious that they differ by the substitution T <-> F and changing 
signs of all variables. Therefore, we get the same Hirota equation for F. We thus conclude that any solution 
to the linear problems (3.10), where the T-function obeys the Hirota equation, provides an auto-Backlund 
transformation, i.e., a transformation that sends a solution of the nonlinear integrable equation to another 
solution of the same equation. In what follows we call them simply Backlund transformations (BT) and 
distinguish two types of them. 

Let us rewrite the linear problems (3.10) changing the order of the terms and shifting the variables: 

T(a+1, s,u)F(a, a, u+1) - T(a, s,u+l)F(a+l, s,u) = T(a+l,s-l,u+l)F(a,s+l,u), 
T(a,s+l,u+l)F(a,s,u)-T(a,s,u)F(a,s+l,u+l) = T(a + 1, s,u + l)F(a-l, s + l,u). (3.16) 

These equations are graphically represented in Fig. 11 in the (a, a) plane. Given polynomials T(a,s,u) 
obeying the Hirota equation, we are going to seek for polynomial solutions for F. 

It is easy to see that equations (3.16) are not compatible if one imposes the boundary conditions for 
F(a, s, u) and T(a, s, u) of the fat hook type with the same K and M. Indeed, applying them in the corner 
point of the interior boundary, one sees that the boundary values must vanish identically. However, it is 
straightforward to verify that equations (3.16) are compatible with the boundary conditions of the following 
two types. The boundary conditions for F(a, s, u) can be cither 

F(a,s,u) = if: (i) a < or (ii) s<0 and a ^ 0, or (iii) a>K-l and s>M, (3.17) 

or 

F(a, s,u) = if : (i) a < or (ii) s < and a ^ 0, or (iii) a > K and a > M + 1 , (3.18) 
which are again of the fat hook type but with the shifts K — > K — 1 or M — ► M + 1. 
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Figure 12: Particular cases of application of the first Backlund transformation (BT1) in the form depicted 
in Fig. 11. 



3.2.1 First Backlund transformation 

The first Backlund transformation (BT1) T — > F is given by the linear problems (3.16) with boundary 
conditions (3.17). Moreover, the linear problems are also compatible with the more specific form of the 
boundary conditions (3.4), where K is replaced by K — 1 and M remains the same. To see this, one should 
check all cases when one of the three terms in the linear equations vanishes, which imposes certain constraints 
on the boundary functions. The first equation at s = (the "brick" in position 8 in Fig. 12), with T(a, 0, u) 
as in (3.2), states that F(a, 0, u) depends only on u + a. Therefore, we can set F(a, 0, u) = Qk-i,m(u + a), 
which at this step is just the notation. Similarly, the second equation at a = (position 5) implies that 
F(0, s, u) depends only on u — s. At s = it must equal Qk-i,m{u), therefore, F(0, s, u) = Qk-i,ai(u — s). 

The consistency of the last two boundary conditions in (3.2) follows from a similar analysis at the interior 
boundaries. Indeed, at a = K — 1, s > M (position 1) the first equation becomes 

Qk,o( u + s + K)Q .m{u -s- K)F(K - 1, s, u + 1) 

(3.19) 

= Qk.o(u + s + K)Q ,m(u - s- K + 2)F(K - 1, s + 1, u) . 

This equation fixes the (u — s)-dcpcndent factor in F(K — 1, s, u) to be Qo.m(u — s — (K — 1)) while no 
restrictions on the (u + s)-dependent factor emerge. One is free to call it Qk-i,o(u + s + K — 1), so that 

F(K - 1, s, u) = Qk-i,o( u + s + k - 1 )Qo-m(u -s- K + 1), M < s < oo, (3.20) 

i.e., the boundary condition for F(a,s 1 u) on the half-line a = K — 1, M < s < oo takes the same form as 
the 3-rd one from eq.(3.4) for T(a,s,u) on the half-line a = K, M < s < oo (see Fig. 10). Finally, one can 
check that the first equation in (3.16) at s — M, a > K (position 6) is consistent with 

F(a, M, u) = {-l) M{a - K+1) Q K -i.o{u + a + M)Q 0>u {u -a- M), K - 1 < a < oo, (3.21) 

and does not bring any new constraints. The sign factor is not fixed by this argument. To fix it, one should 
require consistency with the second Backlund transformation which we consider in the next subsection. 



3.2.2 Second Backlund transformation 

As is seen from eq. (3.18), the transformations generated by the linear problems (3.16) are not able to 
move the vertical part of the interior boundary from the right to the left. For our purpose, we need a 
transformation which would be able to decrease M . Using the duality between T and F explained in the 
end of Sect. 3.2, one can introduce Backlund transformations of the required type. The second Backlund 
transformation (BT2) T — > F* is obtained from eq.(3.16) by exchanging T — > F* and F — ► T: 

F*(a+l,s,u)T(a,s,u+l)- F*(a,s,u+l)T(a+l,s,u) = F*(a+l,s-l,u+l)T(a,s+l,u), 
F*(a,s+l,u+l)T(a,s,u)-F*(a,s,u)T{a,s+l,u+l) = F*(a+l,s,u+l)T(a-l,s+l,u). (3.22) 
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Figure 13: Two equations of the second pair of linear problems in the (s, a) plane. 



These equations are represented graphically in Fig. 13 in the (s,a) plane. As is argued in Sect. 3.2, their 
compatibility condition is the Hirota equation (3.1) for T. If it holds, then any solution F* obeys the same 
Hirota equation and thus provides an auto-Backlund transformation. 

Given a solution T(a, s,u) with the boundary conditions (3.4), equations (3.22) arc compatible with the 
following boundary condition for F*(a, s, u): 

F*(a, s,u) = if : (i) a < or (ii) s < and a ^ 0, or (iii) a > K and s > M - 1, (3.23) 

which directly follow from (3.18) and differ from those for T (3.2) by the shift M — ► M — 1. Moreover, they 
are also compatible with the more specific form of the boundary conditions (3.4), where M is replaced by 
M — 1 and K remains the same. In complete analogy with BT1, one verifies this by applying BT2 on the 
exterior boundaries (positions 8 and 5 in Fig. 12) and on the interior ones (positions 7 and 2). In particular, 
on the half-line s = M — 1, K < a < oo we have 

F*(a,M-l,u) = {-l) (M - 1)(a ~ K) Q K ,o( u + ° + M -1)Qo,m-i{u- a- M +1), K < a < oo, (3.24) 

the sign being uniquely fixed from the second equation of (3.22) applied in position 2. This means that the 
boundary condition for F*(a, s, u) on the half-line s = M — 1, K < a < oo takes the same form as the one 
for T(a, s, u) on the half-line s = M, K < a < oo (see the 4-th equation in (3.4) and Fig. 14). 

Our final goal is to "undress" , using the transformations BT1 and BT2, the original problem by collapsing 
the region where the Hirota equation acts: K — > 0, M — > 0. 



3.3 Hierarchy of Hirota equations and linear problems 

We see now that by applying BT1 or BT2 to a solution to the Hirota equation with the boundary conditions 
(3.4) we can shift the interior boundaries as K — > K — 1 or M — ► M — 1 (Fig. 10 and Fig. 14, respectively). 
In this way we come to the problem for F (a, s, u) or F*(a, s, u) of the same kind as the original problem for 
T(a, s, u). Repeating these steps several times, we arrive at the hierarchy of the functions Tfe iTO (a, s, u) such 
that 

Tk,m(o>, s, u) = T(a, s, u), 

T K _ 1>M (a,s,u) = F(a,s,u), (3.25) 
T K)M _i(a,s,u) = F*(a,s,u), etc., 
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M-1 M S 

Figure 14: Boundary conditions for the second Bdcklund transformation in the (s, a) plane. 



all of them satisfying the Hirota equation 
T fejTO (a, s, u+ l)T fe , m (a, s,u — l)- T k>rn (a,s + 1, u)T k , m {a, s-l,u) = T k<m {a + 1, s, u)T k<m (a -l,s,u), (3.26) 
where k = 0, . . . K and m = 0, . . . , M. The boundary conditions are as follows 

Tk, m (a, s, u) = if: (i) a < or (ii) s<0 and a^O, or (iii) a>fc and s > m, (3.27) 

and 

Tfe, m (0, s,u) = Qk.m(u-s), -oo < s < oo, 

Tfc, m (a,0,u) = Qfe,m(w + a), < a < oo, 

T k!m (k,s,u) = Qk,o(u+s+k)Qo, m (u—a — k), m < s < oo, 

T fc , TO (o,m,u) = (-i)™( Q - fc ) Q fc)0 ( u +a+m)Qo,m(«-a-m), fc < a < oo, (3.28) 

(see Fig. 9). The function Qk,m( u ) ^ s a nxG d polynomial of degree N which depends on the choice of 
the quantum (super) spin chain with the gl(K\M) symmetry. We set <2o,o(w) = 1 (it corresponds to an 
empty chain). The other polynomial functions Q k ,m(u) will be determined from polynomial solutions to the 
hierarchy of Hirota equations. 

The linear problem (3.16) generates a chain of Backlund transformations BT1: 

T k , m {a + l,s,u)T fe _ liTO (a, s,u+ 1) - T k>m (a, s, u + l)T k _ 1>m (a + l,s,u) 

= T ktTn (a + l,s- l,u + l)T fc _i >m (o,s + l,u), 
T ktm (a,s+ L ti + l)T fc _i iTO (a,s,u) - T k<m (a, s,u)T k ^ 1<m (a, s + l,u+l) 

= T k , m (a + l,s,u+ l)T k _ 1<m (a - 1, s + l,u), 

(k = 1, . . . , K). The linear problem (3.22) generates a chain of Backlund transformations BT2: 
Tk, m -i{a + 1, s,u)T k , m (a, s,u + l)- T feim _i(a, s,u + l)T fc)m (a + 1, s,u) 
= Tfc. m _i(a + l,s - l,u+ l)T k>m (a, s + l,u), 

(3.30) 

Tk, m -i(a, s + l,u + l)T fe>m (a, s,u) - T feim _i(o, s,u)T k . m (a, s + l,u+ 1) 
= r fcim _!(a+ l,s,u+ l)T fc , m (a - l,s + l,u) 
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(to = 1, . . . ,M). Note that (3.30) differs from (3.29) only by the "direction" of the Backhand flow: it is 
k — ► k — 1 in (3.29) and to — > to + 1 in (3.30). In the 5D linear space with coordinates a, s, u, k, m, the four 
equations (3.29), (3.30) act in the hyper-planes 



{to 


= const} n {u + s 


+ a 


= const}, 


{to 


= const} n {u — s 


— a 


= const}, 


{k 


= const} fl {u + s 


f a 


= const}, 


{k 


= const} fl {u — s 


- a 


= const}, 



respectively. Thus each of them is actually a dynamical equation in three variables rather than five. It is 
easy to see that all of them can be transformed to the standard form of the Hirota equation (B2) by linear 
changes of variables. 

The two Backlund transformations can be unified in a matrix equation of the type (3.13). Let fk,m (a, s, u) 
be the antisymmetric matrix from the left hand side in (3.13) at the level A;, to: 

/ T k>m (a,s,u-1) -T k>m (a,s+\,u) T ktm (a+l,s,u) \ 

-T k , m (a,s,u-1) T ktW ,(a-l,8,u) T ktm (a,s-l,u) 

T kiTn (a,s+l,u) -T kiTn (a-l,s,u) -T ktfn (a,s,u+1) 

V -T k , m (a+l,s,u) -T k<m {a,s-l,u) T k:Tn (a, s, u+1) / 



Tfc, m (a,s,u) 



(3.31) 



then the Backlund transformations BT1 and (BT2) 1 (the transformation inverse to BT2) in the symmetric 
form are obtained as the first and the second columns of the matrix equation 



T fc!m (a,s,u) 



I Tfc_i )m (a-1, s, u) 
7fc_i, m (a, s + 1, u) 
T fe _i !OT (a, s, u-1) 



7fc,™+i(a-l, s, u) \ 
7fc, m+ i(a, s + 1, u) 
Tk,m+i(a, s, u-1) 



(3.32) 



\ Tfc_i jm (a — 1, s+l,u— 1) Tfe. m+ i(a-l, s+1, u-1) J 



The first and the second equations in (3.29), (3.30) are obtained in the second and the first lines of the 
matrix equation (3.32) respectively. As we have seen above, the rank of the 4x4 matrix T k , m is 2, so it 
has two linearly independent zero eigenvectors. Now we see that they correspond to the two independent 
transformations shifting either k or to. 



3.4 Bilinear equations for the T-functions with different k and m 

In this subsection we derive additional bilinear equations (3.38)-(3.43) for the functions T k>m (a, s, u) in which 
both indices k, m undergo shifts by ±1. They are also of the Hirota type. A special case of them (equation 
(3.45)) is particularly important. It provides a bilinear relation for the Q-functions (the "QQ-relation" ) 
which will be also derived in section 4 by other means. 

All these additional bilinear relations follow from compatibility of equations (3.29), (3.30) with shifts 
in k and in to. To see this more explicitly, we note that these equations admit a remarkable alternative 
representation. Namely, it is straightforward to verify that the first equations in (3.29), (3.30) can be 
identically rewritten as 

Tfc_i, m (a+l,s,n) 
T k -i, m (a, s+l,u) 



T k -i 



■ i,m(a, s, u+l)T k , m (a, s + l,u) _ s 
i, m (a, s + l,u)T k m (a, s, u+1) 



.,-d, 



T k , m (a + l 7 s,u) 
T k , m (a, s + l,u) 



(3.33) 
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Tk,m+i(a + l,s,u) 
T ktm+ x(a, s + l, it) 

while the second equations as 

Tk+i, m (a,s+l,u) _ 
T k+1;m (a+l,s,u) 

T k , m -i(a,s+l,u) _ 
T k>m -x{a+l,s,u) 



T k ^ rn+1 (a, s, M + l)T fc , m (g, s + l, u) ^ du _ ds 
T k ^ m+1 (a, s + l, u)T ktm {a, s, u+1) 



T k , m (a+1, s,u) 
T kim (a,s+l,u) 



Tk+i,m(a, s,u-l)T k . m (a+l, s, u) & -d u -d a 
T k+ i tm (a+l, s, u)T fcjTO (a, s, u- 1) 

Tk, m -i(a, s, u- l)T A . ! „ I (a + l, a, u) e -d u -d a 
7fc im _i(a+l,s,u)T fcjm (a,s,u-l) 



T fc , m (a, s + l,u) 
T feim (o+l, s,u) ' 

Tfc, m (a+l,s,u) ' 



(3.34) 

(3.35) 
(3.36) 



In this form they appear as linear problems for the difference operators in the brackets together with par- 
ticular solutions. One can notice that they again look like auxiliary linear problems for the Hirota equation 
but for a different choice of the variables. For example, equations (3.33), (3.34) should be compared with 
eq. (Bl) from Appendix B, where we identify r with T and choose the variables as p\ = —k, P2 = m, 
P3 = l(u — s) (note that the variables a and i(it + s) are kept constant in the T-functions entering the 
difference operator in (3.33), (3.34)). We see that the two equations coincide with the two corresponding 
linear problems from (Bl) (with Ai = — A2 = 1), with 



Tfc ]m (a+l,s,M) 
Tk, m (a, s + l,u) 



being their common solution. Because they hold for any k,m, the function ?/'fc-i.m+i(a, s, u) can be repre- 
sented in two different ways as follows: 



Jfc-i, m+ i(a, s, M+l)Tfc_i !m (a, s + l, u) 
.Tfc_i jm+ i(a, s+l, u)r fc _i i7n (o, s, u+1) 



„d u -d s 



Tk-i. m (a, s,u+l)T k , m (a, s+l,u) 



_T fc _i, m (a, s + l,u)T k ,m(a,s,u + l) 



Tk-i, m +i(a, s, u+l)Tk, m +i{a, s + l,u) _ e9u _ di 
Tk-i,m+i(a, s+l,u)Tk tm +i(a, s, u+1) 

Opening the brackets, one finds that the terms multiplied by ip ktm (a, s,u) and ipk,m(a, & — 2, u + 2) cancel 
automatically while the terms proportional to ip ktm (a, s— 1, u+1) give a non-trivial relation connecting the 
T-functions with different k and m. It has the form 

T k . m {a, s + l,u)T k+ i, m+1 (a, 8, u+1) - Tk, m {a,s,u+l)Tk+\, m +i{a,s+l,u) 

(3.37) 

= fk,m{a, u+s) Tfe + i jTO (a, s, u+l)T k>m+ i(a, s + l, it) , 

where f k , m {a,u + s) is an arbitrary function of k,m and a, u + s. Comparing it with a similar equation 
obtained in the same way from the other pair of linear problems (3.35), (3.36), one can see that it actually 
depends on the combination u + s — a as well as on k,m (see Appendix C for details). To fix it, we take 
s = in. so the first term of the equation vanishes (at a > k + 1), and use boundary conditions (3.28). This 
fixes the function f k , m to be 1. 

Clearly, eq. (3.37) (with /fc, m (a, u + s) = 1) is the Hirota equation of the form (B2) in the variables —k, 
m and ^(it — s). The Hirota equation implies the compatibility of the linear problems, i.e., the discrete zero 
curvature condition for the difference operators in (3.33), (3.34) holds true. In our situation, it appears to 
be equivalent to the "weak form" of this condition, i.e., with the operators being applied to a particular 
solution of the linear problems. In other words, the compatibility of the linear problems (which in general 
means existence of a continuous family of common solutions) follows, in our case, from the existence of just 
one common solution (cf. [44]). 

There are other equations of the same type. It is convenient to write them all as the following chain of 
equalities: 



T k , m +i(a, s, u+l)Tfc !m (a, s + l, it) 
Tk, m +x{a, s+l,u)Tk, m {a, s, u+1) 



e 



du-d a 
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Tfe,m+i(a, 


s+ l,u)T k+ i tm (a,s,u-\ 


-1) 








l,s, 


u)Tk+i, m +i(a, s 


, u + 1) - T fejTO (a, s,u + 


l)Tfc+l.m- 


-i(a - 1, s, 


a) 








- l,s,u)Tk+i, m (a,s,u-\ 


-1) 






T k ,, n (a — 


1,3, 


u)T k+ i, m+ i(a, s 


- l,u) - T k>m (a, s-1, 




_i(a - 1, s, 


u) 



2fe, m (a, s+ l,u)T fe+ i, m+ i(a, s,u + 1) - T kym (a,s,u + l)T k+1 . m+1 (a, s + l,u) ^ 

(3.39) 
(3.40) 

Jfe,m+i(a - 1, s,u)T k+liTn {a, s-l,u) 

These equations have the same structure. In each equation, one of the variables a, s, u enters as a parameter. 
More precisely, they act in the hyper-planes 

{a = const} n {u + s = const}, 
{s = const} fl {u — a = const}, 
{u = const} n {a + s = const}, 

respectively. The first equality in this chain is already proved. The proof of the other two is straightforward: 
one should pass to common denominator, to group together similar terms and to use equations (3.32). In 
fact this chain can be continued by three more equations of a similar but different structure: 

l _ T k , m (a, s - l,u)T fc+ i, Tra+ i(a, s,u + 1) - T fciTra (a, s,u+ l)r fc+ i ]m+ i(a, s-l,u) 
Tk, m +i(a - l,s,u)Tk+i, m (a + 1, s - 1, u + 1) 



Tk, m {a, s,u + l)Tfc + i, m +i(a + l,s,u) - T k , m (a + 1, s,u)T k+1 ^ m+1 (a, s,u+l) 
Tk, m +i(a, s + 1, u)T k+ i tm {a + 1, s - l,u + 1) 

Tk, m (a, s - l,u)T k +i t m+i(a + 1, s,u) - T k , m {a + h s, u)T k+ i, m +i(a, s-l,u) 
Tk, m +i(a,s,u- l)T k+ i tm (a+ l,s - l,u+ 1) 



(3.42) 
(3.43) 



They are proved in the same manner. There is no need to prove the first equality separately since one can 
just continue the chain of equations (3.38)-(3.40) proving that (3.40) is equal to (3.41). Other forms of these 
equations and more details can be found in Appendix C. Equations (3.41)-(3.43) act in the hyper-planes 

{a — k + m = const} fl {u — s — a = const}, 
{s + k — m = const} n {u + s + a = const}, 
{u — k + m = const} l~l {u + s — a = const}, 

respectively. Therefore, equations (3.38)-(3.40) and (3.41)-(3.43) are actually equations in three variables 
rather than five. They can be transformed to the standard form of the Hirota equation (B2) by linear changes 
of variables. 

We note that equations (3.38)-(3.40) can be written in the following concise form: 

ndv „d a \ ^Wl,m+1 _ ( „d v Tk+l,m \ ( a T k,m+l 



dp ed ±J ^ S±L = edp± ^ ^±^S±L\ . (3.44) 

J- km \ J- km / \ J- km / 

Here T kiTn = Tk tm (a, u, s) and (p, q) stands for any one of the pairs (it, s), (u, —a) and (a, s). 



24 



Restricting the bilinear equations to the boundaries of the fat hook domain, one obtains new equations 
which include the Q-functions. For example, setting a = in eq. (3.37) (or eq. (3.38)) and using (3.28), we 
get the QQ-rclation mentioned in the Introduction: 

Qk,m(u)Qk+l,rn+l{u + 2) - Qk, m (u + 2)Q k+1>ni+ i(u) = Qk,rn+l(u)Qk+l,m(u + 2) . (3.45) 

In the next section, it will be derived by other means. Another new relation, which is of a mixed (TQ and 
QQ) type, is the particular case of eq. (3.43) at a = 0, s = 1: 

7fc+i, m+ i(l, l,u)Qk, m (u) - Tk, m {l, l,u)Qk+i,m+i(u) = Qk,m+\{u - 2)Qfc+i, m (u + 2) . (3.46) 



4 TQ- and QQ-relations 

In this section we partially solve the "undressing" problem for the hierarchy of the T-functions 7fc m (a, s, u) 
and derive the generalized Baxter equations (TQ-rclations) which express T(l,s,u), T(a, 1,m) through the 
Baxter functions Qk,m{u). This is done by constructing an operator generating series for the T- functions 
and factorizing it into an ordered product of first order difference operators, with coefficients being ratios of 
the Q-functions 3 . These operators obey a discrete zero curvature condition which leads to a bilinear relation 
for the functions Qk.m with different values of k, m (the QQ-relation). 

4.1 Operator generating series and generalized Baxter relations 

We start by introducing the following difference operators of infinite order: 

W k<m (u) = > e 

Qk,m(U) 



s=0 



.(«) = E(-l) Thf^lfZ^l (4.1) 

Qk,m{U + 2) 



a=0 



which represent operator generating series for the transfer matrices corresponding to one-row or one-column 
Young diagrams. The denominators are introduced for the proper normalization. Let us show that the 
difference operators 

jj^ ,s _ Qk+l,m{u) Qk,m( U + 2 ) _ & 2dn 

' m Qk+l,m{u+2)Qk,m(u) 

t> („a QkM u ) Qfc,m+i(M+2) 2du 

Vk,m{u) = — — -— — p-r - e . (4.2) 

shift the level indices k, m of the Wk, m (u) and Wk, m (u). Namely, we are going to prove the following operator 
relations: 

W k , m (u) = Uk,rn(u)W k+1 ^ m (u), 
W k ,m+l(u) = Vk,m{u)W k , m {u), (4.3) 

W k +l,m(u) = W k ,m(u)U ktm (u), 

W Km {u) = W ktm¥L (u)V k , m (u). (4.4) 



3 For the bosonic case this was done in [9, 37, 38, 39, 16]. For the supersymmetric case such equations were conjectured in 
[12] (see also [45]). 
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k+1,m 



Ik 



k,m 




k+1,m+1 



▼ Uk,m+1 



k,m+1 



Figure 15: The discrete zero curvature condition Uk,m+iVk+i,- 



To prove the first equation in (4.3), we write 

Wfe-l(u) + e 29 "W k (u) = W k -i{u) + W k (u + 2)e 2a " 



Tn(l,0,« + 1) 



!k-l 



(u) 



E 



T fc _ 1 (l,s, u + S + 1) T fc (l,«-l,u+a+2) 



Qk-i(u) 



Qk{u + 2) 



(4.5) 



(Here and below in the proof we omit the second index m since it is the same everywhere.) To transform 
the expression in the square brackets, we use the first equation of the BT1 at a = (position 3 in Fig. 12): 



T fe (l, s,u)Qk-i(u—s+l) - Tfc_i(l, s,u)Qk{u-s+l) = T k (l, s-l,u+l)Q k -i(u-s-l). 

Shifting it u — > u + s + 1 and dividing both sides by Qk(u + 2)<3/ c _ 1 (it), we obtain 

T fc (l,a,u + s + l)Qfe_i(u + 2) _ T k - 1 (l,s,u + s + 1) = T k (l, s - 1, u + a + 2) 
Qfc(u + 2)Q fe _ 1 (u) Qfc-i(u) Qfc(w + 2) 

Using this, we rewrite the term in the square brackets under the sum in eq.(4.5) as 

Qk(u)Q k -i(u+2) T k (l )S) u + s + l) 



(4.6) 



and continue the equality (4.5): 
W k -i(u) + e 29 -W k (u) 



Qk(u+2)Q k -i(u) 



r Jb _ 1 (i,o,u + i) 



Qk{u) 



)k{u) Q k -i(u + 2) r fc (l,s,u+s + l) e 2 s0u 



Q fc (w + 2)Q fc _ 1 ( U ) 



Qk(«) 



Qk(it) Qfe_i(u + 2) y> r fc (l, s,u + a + 1) e2sau 



) fc (u+2)Q fc _i(u) ^ 



Qk(u) 



(4.7) 



In the last step we have noticed that the s = term of the sum multiplied by the ratio of Q's is just equal 
to Tfe_i(l,0,u+ l)/Qk-i(u). The sum in the r.h.s. is Wk(u), so the first equality in (4.3) is proved. The 
proof of the three other equations in (4.3-4.4) is completely similar. 

Combining equations (4. 3), (4. 4), we see that 

Wk-i,m(u)W k -i, m {u) = Wk, m {u)W k , m {u) = Wk,m+i(u)W k , m +i(u) , 

i.e., the operator Wk, m {u)Wk, m {u) does not depend on k,m. Note that Wo,o(u) = Wo,o(u) — 1 as operators, 
since all the terms in (4.1) are zero except the first one, which is 1 thanks to the "boundary conditions" 
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Q(u+2) 



Q(u) 



Q(u+2) 



Q(u) 



Q(u+2) 



Q(u) 



m 



Figure 16: The Hirota equation for the Baxter functions Qk,m{u) in the (k,m) plane. 

Tofi(l, 0,u+ 1) = Qo,o{u + 2) = 1, Xb,o(0, 1, u + 1) = Qo,o(u) = 1. Therefore, we conclude that the operators 
Wk.m and Wk, m are mutually inverse 4 : 

W ktm {u)W k , m (u) = 1. (4.8) 
In addition, applying equations (4.3), (4.4) many times, we arrive at the following operator relations: 

W kt m = Uj7_ lm ...U^ m Vo, m -i...Vo,o, 

Wk,m = Vq.q ■ ■ ■ ^oTm-l ^0,m ■ ■ ■ £4-1, m, 



(4.9) 



where we have skipped u since it is the same everywhere. Taking equations (4.1), (4.9) at k = K, to = M, 
we obtain the "non-commutative generating functions" for the transfer matrices in the basic representations 
(s 1 ) or (1 Q ): 

ETk ,A/(1, S, U + S + 1) 2sd y_-y - ■ j - - . . 

— ^-7^ n e °" = t/ /f _ 1M (u)...[/ OM (w)Vb,A/-i(w)---Vb,o(w), 

j^ { _ ir T K Ma l,u + a+l) e2adu = ^ o -i (u) ...^_ i(u) ^ m(u) ..^ k _ 1m(u) . (4 . 10) 
^ Qk,m(u + Za + I) 

Expanding the right hand sides in powers of e 2(9u and comparing the coefficients, one obtains a set of 
generalized Baxter relations between T's and Q's. In principle, these formulas solve our original problem: 
they give solutions to the Hirota equation in terms of the Q-functions representing the boundary conditions 
at each level fc, to. 



4.2 Zero curvature condition and QQ-relation 

The Q- functions are polynomials whose roots obey the Bcthe equations. Contrary to the case of bosonic 
gl{K) algebras, the Bethe equations for superalgebras admit many different forms. They correspond to all 
possible "undressing paths" in the (fc, m) plane. Their equivalence can be established by means of certain 
"duality transformations" [14]. 

Here we suggest an easy transparent argument to derive all these systems of Bethe equations and the 
corresponding duality transformations. Namely, we are going to show that the functions Qk,m{u) obey their 
own Hirota equation. Given an undressing path, it immediately produces the chain of Bethe equations. The 
duality transformation is nothing else than the discrete zero curvature condition for the operators (4.2) on 
the fc, to lattice. 

Equations (4.3) imply Wk-i, m +i = ^fc-i,m+iVfc, TO Wfc,m = H-i,i»^-i,m»M which gives the discrete 
zero curvature condition 

Uk,m+lVk+l, m = Vk, m Uk,m- (4-11) 



'Using the determinant representation (2.12), it is not difficult to derive this fact directly from their definitions (4.1). 
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We remark that it looks a bit non-symmetric because V shifts to by +1 while U shifts k by —1. Being 
written through U^ 1 and V, the zero curvature condition acquires the standard symmetric form 

C^+Am = V k+ l, m Uk,l- ( 4 - 12 ) 

As a consequence of it, the following bilinear relation for the Q's is valid: 

Qk,m{u)Qk+l,m+\{u + 2) - Qk+l,m+l(u)Qk,m{u + 2) = Qfc, m +l(u)Qfc+l, m (u + 2). (4.13) 

It was already derived in section 3.4 as a particular case of more general "TT-rclations" (3.38), (3.39) (see 
(3.45)). This is the Hirota equation in "chiral" variables (see Appendix B). Strictly speaking, the zero 
curvature condition (4.11) implies eq.(4.13) up to an additional factor in the r.h.s. depending on k,m 
which remains unfixed by this argument. However, this factor can always be eliminated by an appropriate 
normalization of Qfc,m(w)'s: 

Jk, m 

Qk, m {u) = A k , m J] (u - uf' m) ) , (4.14) 
i=i 

i.e., by choosing the coefficients Ak lTn - Moreover, the result of section 3.4 shows that the boundary conditions 
(3.28) already imply the normalization in which the QQ-relation has the form (4.13) (see the argument right 
after eq. (3.37)). 

The zero curvature condition allows us to represent equations (4.10) in a more general form. Given an 
arbitrary zigzag path ^k,m from (K, M) to (0, 0), the r.h.s. of these equations becomes the ordered product 
of the shift operators along it (see Fig. 17): 

n ( w ) ' 

(x,n)e7ic,A/ 

_ (4.15) 

n 

(x,n)e7/f,j« 

Here we use the natural notation: x is the vector on the lattice with coordinates (fc, to) (k is the vertical 
coordinate and m is the horizontal coordinate!), n = (—1, 0) or (0, —1) is the unit vector looking along the 
next step of the path. In other words, (x, n) is the (oriented) edge of the path from (K, M) to (0, 0) starting 
at the point x and looking in the direction n. In the first equation, the shift operators are ordered from 
the last edge of the path (ending at the origin) to the first one while in the second equation the order is 
opposite. In fact the zero curvature condition implies that equations (4.15) remain true for any path leading 
from (K, M) to the origin provided the shift operators are chosen as follows: 





,(«) 


= U^ hm (u), 


n = 


(-1,0) 


V( x ,„; 


.(«) 


= U k . m (u), 


n = 


(1,0), 


V( X ,„; 


.(«) 


= V^(u), 


n = 


(0,1), 




M 


= Vfc,,„_i(u), 


n = 


(0,-1) 



(4.16) 



Some simple examples of equations (4.15) are given in Section 8. 



s=0 



!K,M 



(u) 



E(-i 

a=0 



a TK,M{a, 1,11 + + 1) 

Qk,m{u + 2o, + 2) 



5 Bet he equations 

The bilinear QQ-relation (4.13) obtained in the previous section (see also section 3.4 for an alternative 
derivation) gives the easiest and the most transparent way to derive different systems of Bethe equations 
and to prove their equaivalence. In a similar way, the generalized TT-relations (3.38)-(3.40) can be used to 
derive a new system of Bcthe-like equations for roots of the polynomials Tfe jm (a, s, it). 
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5.1 Bethe equations for roots of Q's 

To derive the system of equations for zeros of the polynomials Qk.m (4-14), we put u in the Hirota equation 



(4.13) successively equal to Uj' m \ Uj +1 ' m —2, u 
corresponding to a zero of one of the six Q-functions in the equation. After proper shifts of k and to such 
that the arguments of the Q-functions become uf""" 1 ' and uf""" 1 ' ± 2, we get the relations 



(fc,m) 



2, u 



(k+l,rn+l) (fc,m+l) 



and u 



(fe+l,m) 



■2, each 



— Qk,r 



{k,m) 



Q 



k,m U 



(k,m) 



(k, 



2 ) Qk+l,m+l yu) 
2^ Q k -i,m-i (u (k 
2 ) Qk+l,m+i I u 



(k,m] 



!k,m+l [Uj 



<k, m ] 



Qk+l,m (i 



(fe,m) 



(fc,m) 



2 ) Qk-l,m-l \ U 3 



Qk.rn ( 
Qk,m. (u^ 

Qk,m-l ( u j fe ' m) ) Qk+ 
Qk,rn+1 (uf^ 1 ^ Qk-l,m \ Uj 



J 

(k,m] 



Qk,m-1 U 



Qk 



(k,m] 



!k-l,m U 



(k, 



+l,m [ u j 



(k,m] 



= Q 



k—l,m 



(uf'^J Qk,m+1 (u 
( (k,m)\ 



-l,m U 



,{k,m) 
3 

(k.m) 



Qk,m-1 6 

- Qk,m-x ( uf m) + 2] Q k+1 . 



l+1 (uf' m) -2) Qk^m ( 



(k.m' 

4 ' 

(k,m] 



(a) 
(b) 
(0 
(d) 
(e) 
(f)- 



(5.1) 



Here 1 < k < K — 1, 1 < m < M — 1 and j runs from 1 to Jk. m - This is the (over)complete set of Bethe 
equations for our problem. Their consistency is guaranteed by the Hirota equation (4.13). To convert them 
into a more familiar form, let us divide eq. (a) by eq. (c) and eq. (b) by eq. (d). Using also eqs. (e) and (f). 
it is easy to rewrite the system in the form where each group of equations contains the Q-functions at three 
neighboring sites. In this way we obtain the following sets of equations: 



Q 



k—l.m 



(k,m) 



(fc,m) 



(uf' m) 2) Q ktm (uf' m) + 2) Q k+1 . m (uf' m) ) 



lk-X, 



(5.2) 



Qk,m+1 ( u j fc ' Tn) ) Qk,m 



(fc,m) 



2 Q 



k.m — 1 



+ 2 



!k,m+l \ Uj 



(k, 



Qk,m(uf' m) +2) Q fc , m _ X (uf' m) ) 



(5.3) 



>k+l,m \ Uj 



Ik. 



[ uf' m) + 2 



Qk+l,r 



(k, 



2) Qk,m-1 \ U 



(k,m] 



(5.4) 



Qk,m+l (uf^^j Qk-l,m (l 

Q( (k,n 
fc,m+l \ Uj 



(k, 



- 2 



2 Qfc-i, 



/ (k,m)\ 



(5.5) 



These equations are valid at any point of the (fc, to) lattice and do not depend on the choice of the undressing 
zigzag path. The figures show the sites of the (k,m) lattice (k and to are the vertical and horizontal 
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Figure 17: The "undressing" procedure by means of BT1 and BT2 in the (s, a) plane along an arbitrary 
zigzag path Jk,m- This path defines the Cartan matrix which characterizes the basis of generators of the 
Cartan subalgebra. 



coordinates respectively) which are connected by the corresponding Bethe equation. The point (k,m) is the 
one between the other two. The edges of the (k, m) lattice are represented by the arrows which show the 
directions of the transformations BT1 and BT2. For completeness, we also present here two other equations 
derived from (5.1): 



Qk,m+i (uf' m) ) Q Km (uf' m) - 2) Q k+1>m (uf m) + 2) 
Q k , m+ i {uf' m) - 2) Q k , m {uf' m) + 2) Q k+hm (u?- m) ) 




It is clear from the figures that these patterns are forbidden for a zigzag path. 

Let us show how to reduce this 2D array of Bethe equations to a chain. Suppose one fixes a particular 
zigzag path Jk,m from (K,M) to (0,0). Then, for each vertex of the path (except for the first and the last 
ones), one writes the Bethe equations according to the configuration of the path around the vertex. Let us 
enumerate vertices of the path by numbers from to K + M so that the point (k,m) £ Jk,m acquires the 
number k + m. Set 

Jk-\-m 

Qk,m(u) = Qk+m(u) = A k+m JJ (u - uf +m) ) , (fc, m) G 7^,M ■ 

i=i 

Let us also enumerate edges of the path by numbers from 1 to K + M so that the edge joining the (n — l)-th 
and the n-th vertex acquires the number n. (Note that in the course of undressing one passes these edges in 
the inverse order.) To the n-th edge of the path (oriented according to the direction of the undressing, i.e., 
either from the north to the south or from the east to the west) we assign the sign factor according to the 
rule: 

!+l if the n-th edge looks to the south 
— 1 if the n-th edge looks to the west 
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Then the system of the Bethe equations along the path can be written as follows: 

Qn+l + 2>Pn+lj Qn - 2p„+iJ Q n -1 (V^) i +P , lP „ + 1 

Qn+l («f } ) Qn (uf ] + 2 Pn ) Q„_! (u^ - 2p n ) 

where n runs from 1 to K + M — 1. The boundary conditions are Qq(u) = 1, Qk+m(u) = 4>{u). Any chain 
of Bethe equations includes if + M — 1 equations for the roots of K + M — 1 polynomials Qk,m(u) picked 
along a path as in Fig. 17. All the other [(K + 1)(M + 1) - 2] - (if + M - 1) = /Of Q-functions inside the 
K x M rectangle can be expressed through them by iterations of the Hirota equation (4.13). 

This form of the Bethe equations agrees with the general one suggested in [3]. To sec this, let us redefine 
the Q-functions by the shift of the spectral parameter: 



which is equivalent to 



In terms of the roots 



2fc,m(w) = Qk+m (u-k + m) , (5.9) 

Qn(u) = Qn U-f>«) ' ( 5 - 10 ) 



Uj = u) 



4 n) -&- ( 5 - n ) 

of the polynomial Q n (u) the system of Bethe equations (5.8) acquires a concise form 



K +? Qb (uf - K ab ) 

n -)^— -4 = (-1)^ , a=l,...,K + M-l, (5.12) 
lb \ u) + K al 



6=1 W| u i T-rv 6 

where 

A' Q 6 = (pa +Pa+l)5 a ,b ~ - Pa<5a : &+1 (5.13) 

is the Cartan matrix for the simple root system corresponding to the chosen undressing path (see, for 
example, [14]). So it is natural to think of the Kac-Dynkin diagram for the superalgcbras as a zigzag path 
on the (k, m) plane, as shown in Fig. 17. 

Let us give a remark on the duality transformations. There are ^j^jTr different ways to choose the 
undressing path jk,m (Fig. 17), and hence there are as many chains of Bethe equations for gl(K\M) algebra, 
all of them describing the same system but for different choices of the simple roots basis. The transformation 
from one basis to another can be (and sometimes is) called the duality transformation meaning that the 
two descriptions of the model arc equivalent and in a sense dual to each other. For particular low rank 
supcralgebras such transformations were discussed in [46, 49, 50, 51, 52], and for general superalgcbras in 
[14] (see also [45]). In various solid state applications of supersymmetric integrable models (for example, the 
t-J model) this transformation corresponds to the so-called "particle-hole" duality. It is clear that any duality 
transformation can be decomposed into a chain of elementary ones. The elementary duality transformation 
consists in swithing two neighboring orthogonal edges of the path (joining at a "fcrmionic" node of the Kac- 
Dynkin diagram) to another pair of such edges surrounding the same face of the lattice. It corresponds to 
replacing eq.(5.4) at the roots u = u^ k ' m ^ by eq.(5.4) at the roots u = u ^ +1 ' m x ' or vice versa, which induces 
also the subsequent change of the Bethe equations at two neighboring nodes. On the operator level, the 
elementary duality transformation consists in replacing V^~ m Uk,m+i by Uk,mVk+i m m the products (4.10), 
according to the zero curvature condition (4.11). 



5.2 Bethe-like equations for roots of T's 

Actually, roots of all the polynomial T-functions Tfc iTO (a, s, u) obey a system of algebraic equations which 
generalize the Bethe equations (5.8). They can be derived along the same lines using, instead of the QQ- 
rclation (3.45) or (4.13), the bilinear TT-relations (3.38)-(3.43). Fixing an undressing zigzag path jk,m, we 
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set 

Jk + m 

T k , m (a,s,u) = T k+m (a,s,u) = A k+m (a,s) J| (u - uf +m \a, s)) , (k,m) e 7k,m ■ 

3 = 1 

At a = s = the root (0, 0) coincides with uf^ from the previous subsection. Repeating all the steps 
leading to the Bethe equations (5.8), we arrive at the following Bcthe-like equations: 

1+PnPn+l 

(-1) 3 , (5-14) 

(-1)— (5.15) 
(-1) 5 , (5.16) 



T n +i(a, S-Pn+l, W*+Pn+l) Tn(a, 


S+Pn+l, 


u*-p„+i) T n - 


-i(a,s, 




2n+i(a, s,«*) T n (a,s~p n ,u^ 


+Pn) T n 


_i(a, s+p„, u* 


~Pn) 




T n+1 (a+p n+ i,s,u*+p n+ i) T n (a- 


-Pn+1, S, 


W*-Pn+l) r n - 


i(a,s, 


u*) 


T n+1 (a,s,u*) T n (a+p n ,s,u* 


+Pn) T n 


-i(a-p„,s,u* 


~Pn) 




T n+1 (a-p n+1 ,s+p n+1 ,u„) T n (a-\ 


-p„+i,s- 


~Pn+l> M *) ^n- 


i(a, s, 





T n+ i(a,s,«*) T„(a-p„,s+p„,u*) T„_i(a+p„, s-p„, «») 



T„+i(a -p n +i,s, u* 


- Pn+l) In (a, S + p„, 


u* +Pn) T„_i(a+p n , 


s -p„,u») 


r n+ i(a - p n +i,s + p„. 


hi,u*) r n (a+p n+ i,s, 


it* - Pn+i) T„_i(a, s - 


-p„,u* +p„) 


T n+ i(a, s+p n +i,u* 


- Pn+l) T n (a — J3„, s, 


it* +p„) T n _i(a +p n , 


s - p„, u*) 


T n+ i(a — Pn+i, s,u* - 


Pn+i) r„(a+p n+ i,s 


+ Pn+i,u*) r„_i(a, s - 


-p„,it* +p„) 


T„ + i(a, s + p n +i,«* 


- Pn+i) T n (a -p n ,s 


-p„,u*) T„_i(a+p„, 


s,u* +p n ) 



Zn+i(a - Pn+1, 5 +2>n+lj^*) T„(a, 5 - Pn+X,^* ~ Pn+l) T n _i(a + p n , S, ^* + p n ) _ ^_ ^ 1+Pn f n+1 /g 



= (-1)^— , (5.18) 



l+PnP„+l 

(-1— (5.19) 

rt+lj «* — Pn+1 ) Jn^-fc, 5 " Pn, u *J J n-H« Tfti »j «* "T ) 

where u* = uv 1 (a, s). The values of a, s and n are assumed to be such that none of the T-functions is 
identically zero. At the boundaries a = or s = equations (5.14), (5.15) coincide with the Bethe equations 
(5.8). 

6 Algorithm for integration of the Hirota equation 

In this section we develop a general algorithm to solve the Hirota equation (3.26) expressing the functions 
Tk, m {a, s,u) through the boundary functions Qk,m(u) (3.28). We note that it gives an operator realization 
of the combinatorial rules given in [12]. 

6.1 Shift operators 

Our starting point is the alternative representation of the first and second Backlund transformations given 
by equations (3.33)-(3.36) which we rewrite here in a slightly different form. Equations (3.33), (3.34) read 

T k -i, m (a,s,u) = H k - im (a — l,s,u) T kim (a,s,u), (6.1) 

T ktm+1 (a,s,u) = H kim +(a-l,s,u)T k>m (a,s,u), (6.2) 

and equations (3.35), (3.36) read 

Tk+i. m {a,s,u) = H k + tTn (a,s - l,u) T ktm (a,s,u), (6.3) 

T k , m -i(a, s, u) = H ktm -(a,s — l,u)Tk, m (a,s,u). (6.4) 
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H k - m (a,s,u):=— — — — — — — e" % (6.5) 

T kim {a,s,u + 1) T k , m (a,s,u+ 1) 

f> / \ T km+1 (a,s,u + l) T km+1 (a,s + l,u) d g 

£/fc, m +(a,s,w) := — — — — e" % (6.6 

T k>m (a,s,u+ 1) T kim {a,s,u + 1) 

H k +, m (a,s,u) := — -y- + — — e *• " + (6.7 

T k<m {a, s, u - 1) T ktm {a,s,u - 1) 

(«,.,«) - r ^- l( °' B '""i? ) + e"(^^). (6.8) 

T kim {a,s,u-1) T kiTn (a,s,u - 1) 

As we have seen in section 3.4, these equations are equivalent to equations (3.29) and (3.30) which have been 
used to define the Backhand transformations BT1 and BT2. 

The operators introduced above obey, by construction, the "weak" zero curvature conditions 

T k zfi, m ±i = H kTlm ±H k ^ m T k:m = H k ^ m±1 H krn ±T k ^ n , (6.9) 

where H k -, m = H k -. m (a - l,s,u), H k+>m = H k + tm (a,s - l,u), H k>m + = H k>m +(a - l,s,u), H k ^ n - = 
H km -{a,s — l,u) and T k , rn = T ktTn {a, s,u). As is pointed out in section 3.4, they imply the "strong", 
operator form of these conditions: 

#fc=Fi,m± H kT ^ m = H k T tm ±i H k m ± . (6.10) 

The operators H k ± m , H k m ± generalize the shift operators U k _ m , V k<m introduced in section 4. We hold 
the same name for them. Comparing to U k _ m , V k _ m , they act to functions of three variables, not just to 
functions of the spectral parameter u, and involve non-trivial shifts in two independent directions. However, 
the shift operators at a = or s = are effectively one-dimensional since they do not depend on u + s (or 
u — a): 

Hk-, m (0,S,U) 
-fffc,m+(0,S,u) 



H, 



Qk-l,m[U 


— s - 


M) 


Qk-l.m{u - S 


-1) 


Qk,m(u - 


s 4 


1) 


Qk,m{u - S ■+ 


-1) 


Qk,m+l(u - 


- s - 


M) 


Qk.m+l(u - S 


-1) 


Qk,m(u - 


s 4 


1) 


Qk,m(u - S 4 


1) 


Qk+l,m(u ■ 


la- 




Qk+i,m(u 4 a 


+ 1) 


Qk,m(u 4 


a — 




Qk,m{u 4 a- 


-1) 


Qk,m-l( u 


4a 




Qk,m-i(u 4 a 


+ 1) 


Qk,m(u 4 


a — 




Qk.m(u 4 a- 


-1) 



fc+,mV 

(6.11) 

They are functionals of Qk,m(u) only. The first (last) two of them, when restricted to the functions of u — s 
(ii + a), are equivalent to (adjoint) operators Uk, m an d Vfe,m (p km an d V km ) respectively. More precisely, 

#(fc+l)-,m(0>S,u) Q k +l,m(u-S-l) -» U ktm (u-S-1), 



Qk,m(u-S-1) 
1 



Qfc, m +i(«— s-1) 
(-!)"« 



H k>m +(0,S,U) Qk,m(u — S-1) -» Vfc, m (u-S-1), 



Qjfe+i, m (u+a+l) 
(-!)"« 



H k+tm (a,0,u) (-l) a Q k , m {u+a + l) -> Ut m (u+a-l), 



Qk,m{u + a+l) 



H k)(m+1) -(a,0,u) (-l) a Q k>m+1 (u+a+l) -» V£ m (u4a-1), (6.12) 
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where it is implied that the operators in the l.h.s. act on functions of u — s (u + a). 
A simple inspection shows that the shift operators can be written as 

H k -, m {a,s,u) = — — —— - hlJ_ hm {a,s + l,u), (6.13) 

v i \ T k , m+1 (a,s,u + 1) -(i) | -t \ raiA^ 

H k>m +(a,s,u) = — — /i^ +1 (a,s + l,u), (6.14) 

J- k.myO'i s : U -\- 1 ) 
U I \ T k+l,m(a, S,U— 1) -(2) / I -i \ 

H k +, m {a,s,u) = — _ ftfc+i, m (a + M,w), (6.15) 

3- c ^ Tk,m-i(a,s,u-l) *>(2) , 1 ^ 
H ktm -{a,s,u) = — — ? — ftfc ' !(o+ l,s,u) (6.16) 



where 



^fe!L( a ' s ' M ) : = T k , m (a,s,u) (l-e 9u ° 3 ) -, (6.17) 

hf m {a,s,u):= T ktm (a,s,u) (l + e -^ +d °A — -i -. (6.18) 



From this representation it is obvious that they have nontrivial kernels, T ktTn {a, s, u)f k ^(a, u — a + s) and 
Tk,m(a, s, u){~l) a ff, ^(s, u — a + s) respectively, so their common kernel is T kym (a, s, u)(— l) a /j, ^(u — a + ,s), 
where /f*!, are arbitrary functions of their arguments. Modulo these kernels the shift operators H k ± m (a, s, u) 
and H k m ± (a, s, u) can be inverted. We have: 

rr-l / \ rr , , 1 Wi <9 u -d,\-l 1 7fc+i m (a, S, « + 1) 

ff (kfl) -, m (a,-,«) = T fc , TO (a > - + l,«)(l-e« •) Tfcm(a;S + M) r fc , Tra ( W + l) 
^ Tk, m (a, s - j + l,u + j) Tk,m(a, s-j,u + j + 1) 

and 

tj-i ( \ m / i -\ \ / -i I _o„+a„v_i 1 T k m+ i(a, s, u — 1) 

^.(^-(a, -,u)=T fciW (a + !,-,«) (1 + e ->) rfcm(o+M>u) t^K^-D 

= T k , m (a + l,s,u) ± - W^^-l) e"^). (6.20) 

^ Tfe, m (a-j + l, s,u — j) T k , m {a- j, s,u-j —1) 

Equations (6.1) and (6.4) rewritten in the following equivalent form 

T k+1>m (a, s,u) = ^+i)- im (a - Mi") T ktm (a, s,u) , (6-21) 



T ktm+1 (a,s,u) = H k ^ m+l) _(a,s - l,u) T kim (a,s,u) (6.22) 

will be useful for integration of the Hirota equation. 

As it is clear from the explicit expressions (6.19), (6.20), the inverse shift operators acting on the function 
T kl m(a, s, u) in eqs. (6.21), (6.22) are represented by sums of fractions whose numerators and denominators 
are products of T's containing both positive and negative values of the arguments a and/or s. The same 
is also true for products of the shift operators H k ± m (a,s,u) and H km ±(a,s,u) since the operators e~ d " 
and e~ da lower values of a and s. The functions T k-m (a, s,u) are equal to zero at negative integer values 
of s or a according to the boundary conditions (3.27). Therefore, the numerators and denominators of 
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some ratios could simultaneously become zero at some values of s or a. We have to define their values in a 
way consistent with the hierarchy of Hirota equations. One way to do that is to analytically continue the 
T-functions to negative values a = —n + e and/or s = — I + S, where n, I € N and e,(5 S I tend to 0. One 
can straightforwardly verify that the behavior 

T ktm (-n + e,s,u) = 0(e n ), 
T k , n (a,-l + 8, u) = 0(5 l ), 

T k , m (-n + e,-l + 5, u) = 0(e n S l ) (6.23) 

is consistent with the Hirota equation and equations (3.29), (3.30). We then notice that both operators 
-HjT- 1 (a,s,u) and Hrz^_(a,s,u), as well as their products, are nonsingular when acting on the functions 
Tk, m (a, s, u). Actually, the behavior (6.23) is equivalent to the following prescription to define the series of 
the form (6.19) or (6.20) acting on Tfc !m (a, s, u): fractions containing T's at negative values of s and/or a do 
not give any contribution to the sum. Thus, for finite positive values of s and a, the inverse shift operators 
k- m( a ' s ' u ) an< ^ kin- S ' u ) ac thig on T/ S . m (a, s, u) contain a finite sum of nonzero terms. We use this 
prescription in what follows. 

Now, we are ready to present a general algorithm of reconstructing the functions Xfc, m (a, s,u) in terms 
of Qk,m(u) based on equations (6.2), (6.3) and (6.21), (6.22). 

6.2 Integration of the Hirota equation 

Equations (6.3) and (6.21) ((6.2) and (6.22)) are recurrence relations allowing one to express the functions 
Tfe m (a, s, u) in terms of the same functions but with smaller values of k, m and/or a, s: 



T fc , m (a,s,u) = ( ]J Hp-,m( a ~ M,")) T nim (a,s,u), (6.24) 

■p— n+1 
m — 1 

T k , m (a,s,u) = (JJ_H ktq +(a-l,s,u))Tk t i(a,s,u), (6.25) 



q=l 
fe-1 



T^ m (a,s,u) = (JJ^H p+im (a,s — l,u))T ntm (a,s,u), (6.26) 

p—n 

m 

Tk, m {a,s,u) = ( [ H^(a,s - l,u)j Tk,i{a,s,u). (6.27) 

Here k, I, m, n are positive integer numbers such that 

0<Z<m-l, 0<n<fc-l, < k < K, < m < M, 

and UnLz = O m ■ ■ ■ Oi by definition. Substituting n = a, I = into eqs. (6.24), (6.25), and n = 0, I = s 
into eqs. (6.26), (6.27), we obtain: 

K 

T Kfi (a,s,u) = ( H~} Q (a- l,s,u)^ T at0 (a,s,u), (6.28) 

p=a+l 
M-l 

T KM {a,s,u) = Q [ H K>q +(a-l, s,ufj T K ,o(a,s,u). (6.29) 

q=Q 
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(where 0<a<K,0<s< oo) and 



K-l 



T K . M {a,s,u) = ( JJ H p+tM (a, s- l,tt)j T 0iM (a, s,u), (6.30) 



p=0 
M 



T , M (a,s,u) = y _ H Q ^_(a,s- l,u)j T 0tS (a,s,u) (6.31) 
g=s+l 

(where < a < oo, < s < M). This representation allows us to write down the formulas for the whole set 
of the nonzero T-functions which do not belong to the boundaries: 

M-l K 

T K ,Ai{a,s,u) = ( ] [ H K>q + (a - 1, s, «)) ( H~} Q (a - 1, s, u)j T afl (a, s, u) (6.32) 

q— p— a+1 

(where < a < K , < s < oo) and 

K-l M 

T KM {a,s,u) = ( Y[ H p+tM {a,s - l,uj) ( (a, s - 1, u)J T , s (a, s, u) (6.33) 

p— g— s+l 

(where 0<a<oo, 0<s< M). Note that the functions T ai o(a, s, u) and To lS (a, s, u) entering these 
equations are boundary functions: 

T afl {a, s, u) = Q a ,o(u + a + s) (0 < s < oo), T afl {a, s, u) = (-oo < s < 0) (6.34) 

and 

To )S (a,s,u) = (—l) as Qo,s(u— a — s) (0 < a < oo), 2o jS (a, s, u) = (— oo < a < 0), (6.35) 

according to the boundary conditions (3.28). 
Let us put a = 1 in cq. (6.32), 

M-l K 

T KiM (l, s, u) = ( J] (0, s, u)) ( J] i?jT- o (0, s, u)) r ll0 (l, s, u) (0 < s < oo, 1< X), (6.36) 

and s = 1 in cq. (6.33), 

K-l M 

T K>M (a,l,u) = ( Y[ H p+M (a,0,u)) ( [] ^ ( a - °' u )) T QA (a,l,u) (0 < a < oo, KM). (6.37) 

p=0 q=2 

Since the shift operators entering these equations are functionals of the boundary functions Qk. m (u) only, we 
immediately obtain explicit expressions for 7fc jTn (l, s, u) and Tfc !m (a, 1, u) in terms of the boundary functions 
Qk,m{u). 

Let us rewrite the solutions (6.36) and (6.37) in a slightly different form: 

M-l K 

T KM (l,s,u)= ( H H K!q+ (0,s,u)) (Y[H-} o (0,s,u)) T 0fi (l,s,u), (6.38) 

q=0 p=l 
K-l M 

T K , M (a,l,u) = ( [] #p+,mM,")) (II Ho, 1 g -( a >°> u )) ro,o(o,l,«). (6.39) 

Taking into account the explicit form of Xi^o, 

7b o(o,s, u) = 1 at (i) a = or (m) s = and a > 

T 0j o(a, s, u) = at (i) a^0 and s^0 or (ii) s = and a < 0, (6.40) 
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it is easy to see that these solutions can be represented as the generating series 



M-l K 



Y,T K .M(l,n,u-s + n) e "V~-°>) = ( J] H K>q+ {Q,s,u)) (l[H;_\ (0,s, U )), (6.41) 

n—Q q—Q p—1 

oo K-l M 

T k ,m(ti, l,u + a — n) e ~ n{ - du+da) = ( [] H p+M {a, 0, u)) ( JJ (a, O.tt)) . (6.42) 



n=0 p=0 9=1 

Indeed, it is clear that the operator in the r.h.s. of (6.41) is expended in the powers of e 9 "~ 9s as is written 
in the l.h.s., with some coefficients. To fix them, one applies the both sides to the function 2o,o(lj s t u ) an d 
takes into account that e n ^ 8u ~ d ''Too{l, s, u) = To t o(l, s — n, u + n) = unless n = s. The same argument 
works for the second equality. Using (6.12), one can see that the operator relations (6.41), (6.42) are identical 
to (4.10). 

Once the functions Tk m(1> s ? u ) an d Tk m(o>> 1> u ) are constructed, the other T-functions can be expressed 
in terms of Qk,m(u) by either iterating eq. (6.32) with respect to a or eq. (6.33) with respect to s. Therefore, 
setting successively a = 1, 2, . . . (s = 1, 2, . . .) in eq. (6.32) (eq. (6.33))) starting with a = 1 (s = 1), one can 
step by step express Tk,m(ci, s,u) with different a and s in terms of Qk,m{u). Equations (6.32) and (6.33) 
solve the problem of the integration of the Hirota equation for the case of rectangular paths. Using zero 
curvature conditions (6.10), one can easily generalize these equations to the case of an arbitrary zigzag path, 
along the lines of section 4.2. 



7 Higher representations in the quantum space 

In the previous sections, it is implied that zeros of the polynomial function <j>{u) are in general position. 
This corresponds to an inhomogeneous spin chain in the vector representation of gl(K\M) at each site. In 
principle, this includes all other cases such as spins in higher representations in the quantum space. Indeed, 
the higher representations can be constructed by fusing elementary ones according to the fusion procedure 
outlined in Section 2. There, wc have considered fusion in the auxiliary space but for the quantum space 
the construction is basically the same. To get a higher representation at a site of the spin chain, one should 
combine several sites of the chain carrying the vector representations, with the corresponding 6>'s being chosen 
in a specific "string-like" way, and then project onto the higher representation. Before the projection, the 
spin chain looks exactly like the ones dealt with in the previous sections. However, zeros of the function 
4>(u) are no longer in general position. 

The string-type boundary values impose certain requirements on the location of zeros of the polynomials 
T(a,s,u). Indeed, the Hirota equation implies that some of these polynomials must contain similar "string- 
like" factors. From the Hirota equation point of view, the projection onto a higher representation means 
selecting a class of polynomial solutions divisible by factors of this type. In fact, given the boundary values, 
different schemes of extracting such factors are possible. They correspond to different types of fusion in the 
quantum space. 



7.1 Symmetric representations in the quantum space 

To be more specific, consider the simplest case of symmetric tensor representations (one-row Young dia- 
grams) . Fix an integer £ > 1 and consider a combined site i consisting of I sites (labeled by the double index 
as (i, 1), . . . , (i, £)) carrying the vector representation. According to the fusion procedure, the corresponding 
parameters 6> 2 ;, r form a "string": 9 i)r = $i — 2(r — 1), r = 1, . . . ,1. Therefore, the boundary values of the 
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T-functions are given by T(0, s, u) = <ft^ (u — s), T(a, 0, u) = <p^ (u + a), where 

4>j{u) = 4>{u)<j>(u + 2) . . . <j){u + 2{l - 1)) 



1 r( ^ } • 



£ = 0,1,2,. 



(7.1) 



Here </5(it) = YiiLi( u ~ as before. The representation through the T-function is useful for the analytic 
continuation in I. 

A thorough inspection shows that the Hirota equation is consistent with extracting the following poly- 
nomial factor 5 from T(a, s, u) for s = 0, 1, . . . ,£ and a > 1: 



T(a, s, u) 



4>J- S {u + s + a)T(a, s, u) 

N 



e-s)N 



n 



" u — s-\-a- 







L r(3t 



s+a- 



(7.2) 



T(a, s, u) 



Here T(a, s, u) is a polynomial of degree sN (if s = 0, 1, . . . , €). We extend the definition of T(a, s, u) to 
higher values of s by setting T(a,s,u) = T(a,s,u) (s > £, a > 1). Note that the factor in the right hand 
side of (7.2) is a product of functions depending separately on u + s + a and u — s + a. Therefore, if this 
relation between T and T was valid for all values of s, then the T's would obey the same Hirota equation in 
the whole (a, s, u) space (see (2.28)). However, since the definition of T is changed when s > £, the Hirota 
equation breaks down in the plane s = £. It is easy to see that it gets modified as follows: 



T(a, s, u+l)T(a, s, u—1) — [4>(u+s+a — l)] s ' 1 T(a, s+1, u)T(a, s — 1, u) 



4> 



s+(£-s)e(s>^ 



(7.3) 



(u— s) T(a+1, s, u)T(a — 1, s, u) . 



The function 6{s > £) is defined as 6(s > £) = 1 if s > £ and otherwise, so the pre-factor in the right 
hand side equals </>j~(u — s) if < s < £ and 4>^(u — s) if s > £ (at a = 1). The boundary conditions are 
T(o, 0, it) = T(0, s, u) = 1 (a, s > 0). If the point (a, s) belongs to the interior boundary, then the function T 
defined by cq. (7.2) (and by eq. (7.5) below) may contain some additional zeros of a similar string-like type. 



7.2 Antisymmetric representations in the quantum space 

In the case of the antisymmetric fusion (one-column Young diagrams) the parameters 0i :r also form a "string" : 
Q%.r = Si + 2(r — 1), r = 1, . . . ,£. Comparing to the symmetric fusion, this string "looks" to the opposite 
direction, i.e., the sequence of 0's increases rather than decreases. The two types of strings are actually 
equivalent since they are obtained one from the other by an overall shift of <Vs. Our convention here is 
chosen to be consistent with the general case outlined below. 

The boundary values of the T-functions are given by T(0, s,u) = <j>J (u — s), T(a, 0, u) = <f>7 (u + a), where 

(f,-(u) = (f>(u)(t)(u-2)...(f>(u-2(£-l)) 

n pz-u-e, -a (7.4) 
= 2 (WTT V 2 ++J £ = 0,1,2,.... 



5 We note that this factor here and in (7.5) below can be obtained directly from the fusion procedure in the quantum space 
as the product of "trivial zeros" of the fused ij-matrices. 
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A thorough inspection shows that the Hirota equation is consistent with extracting the following polynomial 
factor from T(a, s, u) for a = 0, 1, . . . ,£ and s > 1: 



T(a,s,u) — <t>* (u — s — a)T(a,s,u) 

N 



2 (£-a) 



N 



' u—s—a—b 



n 

i=l 



r 



' u — s-\-a — b 



l + i) 



T(a, s, u) 



(7.5) 



Here T(a, s, u) is a polynomial of degree aN (if a = 0, 1, . . . ,£). We extend the definition of T(a, s, u) to 
higher values of a by setting T(a, s, u) = T(a, s, u) (a > £, s > 1). The modified Hirota equation for T reads 
as follows: 



S.,1 - 



r(a,s,it+l)T(o,s,«-l) - [0 o+ {.g- a )0(a>*)( u + a )J T(o,s+l,u)T(o,s-l,u) 
= [0(u — s — a + l)] 5 ""' f(a+l, s, u)T(a—l, s, u) . 



(7.6) 



The pre-factor in the second term in the left hand side equals 4> a (u + a) if < a < £ and <p e (u + a) if a > ^ 
(at s = 1). The boundary conditions are T(a, 0, it) = T(0, s, m) = 1 (a, s > 0). 



7.3 The general case: remarks and conjectures 

Let us consider the case of a general covariant representation fi at each site of the chain. We construct such 
a site i by fusing \fi\ "elementary" sites with (^-parameters 0i + 2(p — q), where (p, q) G fi. We remind the 
reader that the integer coordinates (p, q) £ Z 2 on a Young diagram are such that the row index p increases 
as one goes from top to bottom and the column index q increases as one goes from left to right, and the top 
left box of n has the coordinates (1,1). 

Given a diagram fj,, we define the polynomial function 

N 

M«)=n n («-2(p-g)-0i)= n <^-2( P - g )), (7.7) 

then the boundary values of the T-functions are: 

T(0,s,u) = (f>v(u — s) , T(a,0,u) = <j)^(u + a) . (7.8) 



Let /Lt(a, s) be the Young diagram obtained as the intersection of /i and the rectangular diagram (s a ) 
(Fig. 18): 

fi(a, s) = fj, R (s a ) . 

For brevity, we will sometimes denote the cut diagram p,(a, s) simply ft dropping the dependence on a, s. 
Let us introduce the polynomial function 

<t>»\ iiM {u)= ^ {U) = [] 0(«-2(p-g)), (7.9) 

(p , g)eMMa , s) 

where the product goes over boxes of the skew diagram fj, \ p,(a, s). If fi is contained in the rectangle (s a ), 
then we set </> jU \ j5 ( q , s ) (u) = 1. 

Now we are ready to present our first conjecture. We expect that the projection onto the representation fi 
in the quantum space means, for the Hirota equation, that we consider the solutions such that the polynomial 
T(a,s,u) is divisible by the polynomial 0^\/j( a ,s) : 

T(a,s,u) = (j>n\ji( atS )(u-s+a)f(a,s,u) . (7.10) 
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s 



M(a,s) 




|j\|j(a,s) 



Figure 18: The intersection of Young diagrams y, and (s a ). 



For one-row or one-column diagrams this formula yields equations (7.2) and (7.5). Presumably, this conjec- 
ture can be proved by means of the technique developed in [36]. If the point (a, s) belongs to the interior 
boundary, then the function T defined by eq. (7.10) may contain some additional zeros of a similar string- like 
type. 

We note that the functions T(a,s,u) defined by (7.10) obey the modified Hirota equation: 

T(a, s, u + 1) T(a, s,u — l) 



(f>(u+s+a-2j-l) 

1=^+1 



T(a, s+1, u)T(a, 8 — 1, u) 



(7.11) 



Here the products like T[ r 



</)(u-s-a+2j+l) 

3=1*0+1 

ire understood to be 1 and 



T(a+1, s, u)T(a — 1, s, u) 



y! s = min (/x' s) a) , ^ a = min(^ a ,s) 

are respectively the lengths of the s-th column and the a-th row of the diagram p, — /2(a, s). The boundary 
conditions are T(a, 0, u) = T(0, s, u) = 1 (a, s > 0). Note that the pre-factors in the modified Hirota equation 
are equal to 1 if the rectangle (s a ) is contained in the diagram [i. In the opposite case, when the diagram \i is 
contained in the rectangle (s a ), the functions T(a, s, u) coincide with T(a, s, u) and the pre-factors are again 
equal to 1. Given equation (7.10), the derivation of the modified Hirota equation for T's is straightforward. 
We present here two simple identities which appear to be useful: 

(f>p,(a, S +l)(u) 



(f>p,{a,s)(u) 

<ftft(q+i,s)( M ) 

<f>p,{a,s){u) 



(u+2s) <j){u+2s-2) ... cj)(u + 2s-2(fi' s+1 -l)) , 



(u-2a)<t>(u-2a+2) ... (f>(u-2a+2(p a+1 -l)) 



(7.12) 



The next challenge is to find out how the polynomial factor extracted from the T-functions behaves 
under the Backhand transformations. Here is our second conjecture. At each step (fc, m) of the chain of the 
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transformations BT1 and BT2 the same relation (7.10) holds, 

T ktm (a,s,u) = <t> Mh , m \ji kim (a,8)(u-s+a)T kim (a,s,v,) , (7-13) 

where the diagram /ik,m is obtained from /x = [Ik.m by cutting off K — k upper rows and M — m left columns. 
If K — k exceeds the number of rows in the diagram /x, or M — m exceeds the number of columns in /x, then 
we set Hk.m = 0- In other words, the transformation BT1 cuts off the upper row while BT2 cuts off the 
left column. The coordinates (p, q) on the diagrams fik,m are such that the top left box of any diagram has 
coordinates (1,1). 

At last, we would like to remark that instead of the function </> m (m) one could use the function 

4>M= II 4>{u + 2{p-q)) (7.14) 

and arrive to similar formulas. This can probably be explained by invoking the representation theory of 
(super) Yangians. It suggests 6 that the function 4>^{u) corresponds to the representation associated with 
the usual Young diagram /x while the function </^(u) comes from the representation associated with the 
"reversed" diagram regarded as a skew diagram, see the first reference in [36] . (The reversed diagram is 
obtained from /x by inversion with respect to the left top corner.) This point needs further clarification. 



7.4 Bethe equations with non-trivial vacuum parts 

It is known that when "spins" in the quantum space belong to a higher representation of the symmetry 
algebra, Bethe equations (5.2)-(5.5) acquire non-trivial right hand sides (sometimes called vacuum parts). 
We are going to show that they actually follow from the equations with trivial vacuum parts if one partially 
fixes the roots of the polynomials Qk.m{u) in a special way. Specifically, we set 

Qk,m(u) = (t>n k , m (u) Qk,m(u) , (7-15) 

which can be regarded as an ansatz suggested by the fusion procedure. Note that it is the particular case 
of eq. (7.13) at a = or s = 0. Accepting this, we are going to substitute it into the QQ-relation to get an 
equation for Q's. This allows us to derive the system of Bethe equations for the roots of Q's. The derivation 
itself does not depend on the validity of the conjectures given above. 

To proceed, we need some more notation. Given a Young diagram /x, let = /x^ be the number of its 
rows or, equivalently, the length of the first row of the transposed diagram /x'. The short hand notation 

lk,m = l{fJ>k,m) . l'k, m = l (^k,m) ( 7 - 16 ) 

is convenient. In other words, lk,m x l'k m IS ^ ne niinimal rectangle (of height lk, m and length l' k m ) containing 
the diagram fXk,m- It is obvious that if Hk, m $j then (sec Fig. 19) 

lk,m = K, m )l = MM-m+1 -K + k 

(7.17) 

l 'k,m = (Mk,m)l = VK~k+l ~ M + m ■ 

Now we are ready to substitute (7.15) into the Hirota equation for Q's (4.13). We have: 

= ^»_i,„(«-2)^ im («) 

= <?Wm-l (« + 2 ) 0/" m («) , 

6 We are grateful to V.Tarasov for a discussion on this point. 
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M-m 



M 



K-k 



K 




Figure 19: Reduction of the Young diagram fj. to fJ, k , m under K —k applications of BT1 and M —m applications 
of BT2. The full diagram \x is depicted by thin lines (including the axis), the reduced one, fJ-k,m> by the thick 
lines. 

provided fJ,k,m is not empty (otherwise the right hand sides are equal to 1), where the "string polynomials" 
4>f(u) are defined in (7.1) and (7.4). Using these obvious identities, it is straightforward to obtain the 
QQ-rclation: 

H U + 2l 'k, m ) Qk-Urn-l(u) Qk, m (u + 2) 



u-2Z fe>m +2 



>k-l,m-l 



(7.18) 



= Qk-l,m(u) Qk,m-l(u + 2) . 

In this form it is valid if fj, k , m 7^ 0, otherwise the functions Q obey the standard Hirota equation (4.13). In 
other words, the Hirota equation gets modified in the region shown in Fig. 20. The boundary of this region 
in the (fc, m) plane is exactly the boundary of the diagram /i. The functions Qo,o and Qk,m are fixed to be 
1: Qo,o(w) = Qk,ai( u ) = !■ 

The Bethe equations are derived from (7.18) in the same way as in section 5. They acquire non-trivial 
right hand sides which can be compactly written in terms of the quantities lk, m and l' k m : 



Qk-i, m (uf ,m) ) Q Km {uf' m) 2) Q k+1 , m {uf - m) +2) <f> (uf'^ + 2l' k>m ) 



ffe-1, 



(uf m) - 2) Q k , m (uf' m) + 2) Q k+1 , m (uf' m) ) (uf m) + 2l' k+1>m ) 



(7.19) 



Q k , m+ i (uf m) ) Q Km (uf' m) - 2) Q k , m -i (uf m) + 2) (uf > m) - 2l k>m+1 ) 



(uf' m) - 2) Q Km (uf' m) + 2) Q k , m -i («f m) ) 4> (t 



(fc,m) 



2/ 



k.rn 



(7.20) 



Qk+l.m ( u f'™J 



h + i,m (uf' m) +2)Q k , m . 1 (uf' m - 



<f)(u 



(fc,m) 



2/ 



<j>(uf' m) -2l k>m 



(7.21) 
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Figure 20: The domain (shadowed region), where the Hirota equation for Q's is modified in the rectangle 
< k < K, < m < M. The diagram (i is reflected with respect to its vertical boundary. 



Q 



(^ m) )Q,-i, m («f' m) -2) _ 4(1 



(fc,m) 



21 



Qk. 



m+1 



( (k,f 



2 I Qk-i, 



uf' m) + 2l' k , m 



(7.22) 



These equations are listed here in the same order as equations (5.2)-(5.5). For empty diagrams, l km and 
l' k m are put equal to 0. Using formulas (7.17), we can represent the right hand sides in a more explicit but 
less compact form: 



Qk-x,m (uf'^) Qk,m (uf' m) - 2) Q k+1 , m (uf' m) +2) <f> (uf - m) +2m-2A/ + 2 MA -_ fc+1 ) 

Q fe -i. m {uf' m) 2) Q k , m (uf m) + 2) Q k+Um (uf> m) ) ' <j>(uf m) +2m-2M + 2vL K - k 



(7.23) 



Qk, m+ x (uf' m) ) Q k , m (uf' m) - 2) Q Km ^ (uf m) +2) <f> (u?> m) -2k+2K-2n' M _ m ) 

Q k , m+ i (uf- m) - 2) Q Km (uf m) + 2) Q kim _ x (uf- m) ) ~ ~J[~ 



u) ■ '-2k+2K-2n' M _ 



M-m+l 



(7.24) 



Q k+ l,m (uf> m) ) Q k , m -1 (uf m) +2) <f> (uf ' m) + 

Q k +l,m ( 



uf' m) + 

> + 2) Q fc , m _r(«? 



</>(wf'" l; +2m-2Af + 2 MA -_A 



(j)[uf' m) -2k+2K-2^ 1 



M-m+l 



(7.25) 



m+l 



(uf m) ) Q fc -i, ro (uf m) - 2) ( 



(uf' m) - 2) Q fc _ 1>m ( u f ( u J fc ' m) +2m-2Jlf +2 MJf _ fc+1 ) 



(7.26) 



These equations are the building blocks to make up the chain of Bet he equations for any undressing path. 
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For example, the chain of the Bethe equations for the simplest path (K,M) — ► (0, M) — ► (0,0) is as 
follows. Moving down from (K, M) to (0, M), we have the equations 



Q k -i,M (uf> M) ) Q k , M - 2) Q k+liM (uf> M) + 2) (uf ^+2^_ feH 



Q k -i,M (4 fe ' M) - 2) Q k<M [uf>^ + 2) Q k+1M (uf^) <$> (uf' M) +2^ k ) 



(7.27) 



where k = 1,2, . . . ,K — 1. They agree with the chain of Bethe equations presented in [47, 48] for the bosonic 
case. At the corner, the equation is 



Qi,M (>? M) ) Qo,M-x + 2) <t> (uf^+2^) 



Finally, moving to the left from (0, M) to (0,0), we have the equations 

Qo, m+ i (uf"°) Qo, m («? ,m) - 2) Qo, m -i (uf- m) + 2) (4°' m) +2^-2^ 



(7.28) 



M—m 



Qo, m+ i (uf' m) - 2) Q 0>m ^ + 2) Qo,™^ ( U f < m) ) («<° - m) + 2A'-2^ 



M-m+l 



(7.29) 



where m = 1, 2, . . . , M — 1. One can see that the differences of arguments of the 0- functions in the numerator 
and the denominator are equal to the doubled Kac-Dynkin labels (A2) corresponding to the diagram fj,. 

One can also write down the chain of equations for an arbitrary path, similarly to eq. (5.8). However, 
the general form of the right hand sides is not very illuminating. 

The last remark concerns the form of the ansatz (7.15). Instead of (7.15) one could use a similar ansatz 

Qk,m(u) = 4>^ m (u+2{K-k)-2(M-m)) Q k , m (u) (7.30) 

with the function </> M (w) defined in (7.14). This leads to a similar system of Bethe equations with non-trivial 
vacuum parts. There is a one-to-one correspondence between solutions of the both systems. 

8 Examples 

8.1 Baxter and Bethe equations for gl(l\l) 

Let us consider in detail the g/(l|l) case. In this case all elements of transfer matrix T\ t x{a, s, u) = T(a, s, u) 
lay on the boundaries of the fat hook domain, as we see in Fig. 21. We read off from this figure: 



r(o,o,«) 




N 

1 = fa) = n (« - OA 

3=1 


T(a, 0, u) 


= Qi,i(u 


+ 0), a>0 


T(0,s,u) 


= Qis(u 


— s), —00 < s < 00 


T(a, l,u) 


= (-1)°- 


1 Qi,o(u + a + l)Q ,i( u - 


T(l,s,u) 


= Qi,o(u 


+ s + l)Q ,i(u-s-l), 



-1), 0>1 

s>l. (8.1) 
Using the Hirota equation (4.13) we get the only nontrivial QQ- relation: 

Qo I o(u)Qi,i(u + 2)-Qi I i(u)Qo,o(u + 2) = Qo,i(«)Oi,o(« + 2), 
which gives, together with Qo,o(u) = 1 ; 

Qi,o(«)Qo,i(« - 2) = 0(«) - 0(u - 2). (8.2) 
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This relation allows us to exclude, for example, Qo,i or <5i,o from eq.(8.1) and to obtain, in particular, the 
Baxter TQ-relations: 

T(l,s,u) = ® oA . U = * ~ ]\ {4>{u + 8+1)- 4>{u + a - 1)) , 

T(l,s,u) = TT^j—^—T^l ( < f ) ( u - s + 1) — <t>{u - s - 1)) . (8.3) 
From the second relation, using regularity of the left hand side, we obtain the following Bcthc equation: 

<t>(u?' 0) ) 

1 = 10 ) (8-4) 

It is a free fermion equation. We can also obtain it from the undressing procedure Qi t i — > Qi.o — * Qo,o = lj 
putting u = u^ 1 ' ) + s — 1 in eq.(8.2). For another undressing procedure, Qi,i — > Qo,i - > Qo.o = 1, or by 
canceling poles at u = u^ ' 1 ^ — s + 1 in the first equation in (8.3), we get the Bethe equation 



'(uf'^+2) 



Note that the Bethe equations (8.4), (8.5) do not depend on the representation parameter s whereas the 
eigenvalues T(l, s, u) do depend on it. In principle, nothing prevents this parameter to be continued analyti- 
cally to the complex plane. This is precisely the continuous label of typical (long) irreducible representations 
[21] of the superalgcbra in auxiliary space. 

Let us consider higher rectangular irreps in the physical space. To cover representations with an arbitrary 
spin £ (one-row diagrams with £ boxes), we take a special form of the function cf>(u), as in cq. (7.1): 

4>{u) = ip(u)ip(u + 2)...<p(u+ 2£)ip (u + 2(1 - 1)) , (8.6) 

describing a "string" of length £. (Comparing to the previous section, we have changed the notation slightly 
denoting the polynomial function with roots 0i by f(u) = Yij( u ~ According to (7.15), we have: 

Qo,i{u) = Qo,i(u) 

Qi,o(«) = p(u)¥>(u + 2)...p(u + 2(*-2))Qi,o(u). 
From the first equation in (8.3) we obtain the transfer matrix Ti^l, s, u) = T(l, s, u): 

T(l, a, u + 1) = ggii^— ij ^(u + s + 2£)- <p(u + a)} TT ip(u + a + 2j). (8.7) 
Qo,i( u + s ) jJi 

The last factor represents trivial zeros of T(a, s, u). It has the same origin as the one in cq. (7.2). However, 
one should note that because T(l, s, u) lives on the interior boundary, this factor contains more zeros than 
the one in eq. (7.2) (see the remark in the end of section 7.1). The Bethe equation (8.5) becomes: 

1 = Li, (8.8) 

<p(uf ' 1} + 21) 

It is easy to generalize these results to the case of an arbitrary physical spin £i at each site I = 1, 2, . . . , N 
of the chain. In this case we take 



4>{u) = ]J{u - 0j)( u _ 0, + 2) . . . (u - 0, + 2£i - 2). 
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Qi,i(u+a) 



Qi,o(u+a+1)Qo,i(u-a-1) (-1) 



a-1 



Qi,o(u+s+ 1 ) Qo,i(u-s- 1 ) 



Ch,i(u-s) s 

Figure 21: Boundary conditions for the transfer matrix T(a, s, u) = ri,i(o, s, it) for the <?^(1|1) superalgebra. 



This yields the transfer matrix 



T(l,s,u + 1) 



Qa,i(u - s) 
Qo,i( u + s ) 



N 



\(u-6i + s + 2li)- \{u 



N \ N ti-1 

+-) niK«- 

1=1 / 1=1 j=0 



+ s + 2j). (8.9) 



The last factor represents the trivial zeros which can be absorbed by normalization. Introducing a new 
rcnormalizcd function T' with the trivial zeroes removed, 



n e,-i 

T(l,s,u+l)=T'{l,s,u+l) ]J Y[(u-6i + s + 2j) 

1=1 d = i 



(8.10) 



we obtain, redefining the parameters 9i = &i + ii, 

N 



rrUi , i\ Qo,l{u-s) 

1 (1, s, u + 1) = — 



Qo,i{u + s) y-f-j^ 



\[u-ei + s + i l 



n ( u -Si+s-it 

1=1 



(8.11) 



which is essentially the same transfer matrix eigenvalue as for the rational limit of the equations (A. 9), 
(A. 10) from [53], with the definition of the gl(l\l) S*-matrix (3.1), (3.2) from [54]. The Bethe equations are 



N (04) n i i 

u) — 6i + s + . 



n 



?; (oa) - 6, + s 



.12) 



For the complete comparison one should exchange in these papers the parameters £ — * s, s — > r of the auxil- 
iary and quantum spaces and take the rational limit. Some simple shifts in the arguments and definitions of 
Bethe roots are also necessary. Note that unlike [53], where they are the soliton charges, both representation 
labels I and s can be now considered as continuous parameters of the typical representation of gl(l\l). This 
corresponds to the limit of large charges and large period of these soliton charges. 



8.2 Baxter and Bethe equations for gl(2\l) 

In the gl(2\l) case 7 , the transfer matrices T2 j i(a, s, u) = T(a, s, u) are on the boundaries of the fat hook except 
T(l, s, u) in the middle row, as we see in Fig. 22. There are six Q-functions Qk,m(u) (k = 0, 1, 2; m = 0, 1), 
two of them being fixed by the boundary conditions. The T- functions are expressed through them in the 

7 The construction of the Baxter Q-operators and TQ-relations for the models based on gl(2\l) (and U q (gl(2\l))) were recently 
discussed in [55], [56]. We thank Z. Tsuboi for bringing these works to our attention. 
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a 



Q2,i(u+a) 



Q2,o(u+a+1)Qo,i(u-a-1) (-1)a-2 
/ 

Q2,o(u+s+2) Qo,i(u-s-2) 




XT T 2 ,i(1,s,u) 



Q2,i(u-s) & 

Figure 22: Boundary conditions for T(a, s, it) = T2,i(a, s, u) for the gl(2\l) superalgebra. 



following way: 



N 



T(0,0,«) 

T(a, 0, u) 
T(0,s,u) 
T(a, l,u) 
T(2,s,u) 



= Q 2il («) = <f>(u) = Y[(u-6j) 
j'=i 

= <22,i(« + a), a > 

= Q2,l{u — s), — OO < S < OO 

= Q 2 ,o(w + a + l)Qo,i("-a-l)(-l)' 



a > 2 



: j0 (w + s + 2)Qo 1 i(m 



2), 



s > 1. 



(8.13) 



8.2.1 Baxter equations for the Kac-Dynkin diagram 1 (g) 

(Hero and below the subscript 1 means the fundamental representation in the quantum space.) 
The operator generating series (4.15) reads 

^ T(ffl,l,M + ffl+l) 2adu _ J J _ 

^ Q 2il ( U + 2a + 2) G ~ U ^ U ^o ~ 



a=0 



M u ) 2d u \ ( Q2,o( u ) Qi,o(w + 2) 2d \ / Qa,o(w) Q2,l( u + 2 ) 



e*«- ) ( _2±p^i - B ^ ,v " - - ) ( „"";" v "' ,^ r M '' - e 2Su ) (8.14) 



,Qi, (u + 2) yVQ2,o(w + 2) Qi, (tt) ;V<32,o(w + 2) Q a ,i(«) 

In particular, 

T(2,l,tt + 3) _ Q 2 ,q(m + 6) <ft(u + 4) Q 1>0 (tt + 4) Q 2 , (n + 6) <j>{u + 2) 

</>(« + 4) " Q 2 , (u + 4) 0(u + 6) Oi,o(it + 2)Q a ,o(« + 4)^(u + 6) 

Qi,o(") Q 2 ,o(m + 6)4>(m + 2) Q 2 , (n + 6) 
Qi,o(u + 2) Q 2 , (u + 2) 0(u + 6) Q 2 ,o(u + 2) cf>(u + 6) ' 

From this equation and equation 

rrtr, \ Q2,o( u + s + 2 ) Wo n L11 
T 2, s, u) = - — —— T 2, 1, u - s + 1 

Q 2 ,o(w - s + 4) 

which follows from cq.(8.13), we obtain 

T(2,s,u) Q 2 ,o(u + s + 2)4>(u-s + 2) Qi, (u - s + 2) Q 2 ,o(u + s + 2) 4>{u - s) 



-i 



(8.15) 



4>(u-s + 2) Q 2fi {u-s + 2) cj)(u-s + 4) Qi,o(u - s) Q 2 ,o(u - s + 2) <f>(u - s + 4) 

Qi,o(m-s-2)Q 2 ,o(u + s + 2) Q 2]O ( u + 5 + 2)0(n-s-2) 

Qi,o(« — s) Q 2 ,o( u - s ) 0(w-s + 4) Q2 : o("- s ) 0(u-s + 4) 



(8.16) 
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Note that the representation index s can be treated as a continuous parameter. It corresponds to the 
continuous label of typical representations of the <?^(2|1) superalgebra. 



8.2.2 Spins in higher irreps in the quantum space 

As an illustrative example to section 5, as well as for the purpose of comparison with the result of [6] (eq.(El) 
obtained from the specific S'-matrix by the algebraic Bethe Ansatz), we consider here spins in the quantum 
space in the irrep (r 2 ) (Young diagrams with two rows of length r). At the end, the number r will be treated 
as a continuous label. 

Let us introduce new definitions, following eq. (7.7) of the previous section: 

= 0a,i(«) = tp(u - 2) [ip 2 {u)ip 2 {u + 2) . . . if 2 (u + 2r - 4)] if (u + 2r - 2) , (8.17) 

<foo (u) = <p(u - 2) [<p 2 (u)<p 2 (u + 2) . . . Lp 2 (u + 2r - 6)] ip (u + 2r - 4) , (8.18) 
0i jO (u) = tp(u)tp{u + 2)...tp{u + 2r-4). (8.19) 
Here <p(u) = Y[j( u ~ Then, according to cq. (7.15), 

Oa I o(«) = 02,o(«)Qa I o(«), (8.20) 

Oi,o(«) = 0i,o(«)Oi,o(«)- (8-21) 
In this notation, eq.(8.16) takes the form 

T'{2,s,u) = QM u + a + 2) 
Q 2 . Q {u -a + 2) 

Qi,o( u — s + 2) Q2,o(u + s + 2) 4*ifi(u — s + 2) <f)(u - s) 
Qi,o(u-s) Q 2fi {u-s + 2) <t>i,o(u-s) <f>(u-s + 2) 

Qi,o(u - s - 2) Q2,o(u + s + 2) <f>\fi(u - s - 2) 2 ,o(" - s + 2) cf>{u - s) 
Qi,o(u-s) Q2,o(u-s) 0i,o(u-s) <fa,o(u-s) <f)(u-s + 2) 

Q2,o( u + s + 2) 2 ,q(m -a + 2) <f>{u - a - 2) 
Q2,o{u-s) 02,o(w-s) (j>{u-s + 2)' 

where we redefined the transfer matrix extracting trivial zeros: 

K ' ' ' V ' ' ; 0(u-s + 4) 2iO (u-.s + 2) 

= T'(2,s,u) [ip{u-s)(p(u-s + 2)<f> 2t0 (u + s + 2)}. (8.22) 

For the reason already discussed in the case of gl(l\l), T" is not equal to T of eq.(7.10). 

Calculating the ratios of 0-functions in eq.(8.22), we find for the following result for the transfer matrix 
in physical irrep (r 2 ) and the auxiliary irrep (s 2 ): 

Tr(2,s,u) 



Q2,o( U + 


s 


+ 2) 








Q2fl( U ~ 


s 


+ 2) 








Qi,o( u - 


s 


+ 2) Q 2 ,o(u + 


s 


+ 2) tp(u — s 


-2) 


Qi,o(" 




s) <?2 : o(u- 


s 


+ 2) <p(u - s - 


f 2r) 


<9i,o(u - 


s 


-2)Q 2 ,o(« + 


s 


+ 2) ip{u — s 


-2) 






S) Q2fi( U 




s) ip{u - s - 


Y2r) 


Q2fl(u + 


S 


+ 2) ip(u- 


s 


- 4) ip(u - 


8-2) 



Q2,o(u-s) <p(u-s + 2r-2)ip(u-s + 2r) 



(8.23) 
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Using the notation of eq.(5.9) Qk+m{ u ) — Q k,m{u + k — m) and setting <p(u) = ip(u + r — 1) wc bring eq.(8.23) 
to the form 



T' r {2,s,u) = 



Qs(u + s) 
Q 2 (u - s) 

Qi(u — s + 1) Q 2 (u + s) (p{u - s - r - 1) 
Qi(u - s - 1) Q 2 (u - s) <p{u -s + r + 1) 

Qx{u - s — 3) Q 2 (u + s) y(« - a - r — 1) 
Qi(u - s - 1) Q 2 (u -s-2)<fi(u-s + r + 1) 

Q2(u + s) <p(u — s — r — 3) ^?(u — s — r — 1) 
Q 2 (u - s - 2) <p{u - s + r - 1) <p(u - s + r + 1) 



.24) 



This transfer matrix eigenvalue coincides, after shifting r — > r — 1 and making an easy generalization to the 
inhomogeneous chain (i.e., to spins in arbitrary irreps at each site of the chain), with eq.(El) taken from [6] 
(see Appendix E, where their result is rewritten in our notation). Note also that the spin label r, as well 
as s, can be treated as continuous parameters here. Hence our method is general enough to describe the 
transfer matrices in all possible typical and atypical irreps in quantum and auxiliary spaces. 

Canceling the poles at u = + s — 2 and u = + s in eq.(8.23), we write the following Bethe 

equations: 



Lp(uf> 0) +2r-2) Qi,o( U f 0) ) 



Qi,o(^- 0) + 2) g 2 ,o(^f 0) ) 

Qi.o («?' 0) - 2) Q 2 , («f 0) + ' 



or, in the notation of eq. (8.24). 

(p (uf + r + l) Qi (4 



1 



(«f - r - l) Qr (uf l) 

Qi + 2) Q 2 (uf 



Qi («?> - 2) Q 2 ( u f + 1) ' (8 ' 26) 

As an example, let us also consider two other equations for different choices of the Kac-Dynkin diagram. 



8.2.3 Bethe equations for the Kac-Dynkin diagram Qi 

Let us express T(l, s, u) through Q's using eq.(4.10): 
^ T(a,l,u + a+l) 2adu _ - _j. j - _ 



h,i(u + 2) g2a \ ( Q 1A (u) Q ,i(m + 2) _ e2a \ ( Q 2 ,i{u) Q 1A (u + 2) _ &2du 



Qo,i(«) / \Qi,i(u + 2) Q ,i( u ) ;V ( 92,i(" + 2) Qi,i(«) 



.27) 
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In particular, 

T(l,l,u + 2) Qi,x{u) Q ,i( u ) <f>(u + 2)Q 1A (u + 4) Q 0A (u) cf>(u + 2) 



(8.28) 



<j>(u + 4) Q M (u + 2) Q 0A {u + 2)(f>{u + 4)Qi t i(u + 2) Q 0)1 (u + 2) 0(u + 4) ' 

Canceling the poles in eq.(8.28), it is straightforward to write the following Bcthc equations for gl(2\l) 
superalgebra: 

0( U ; M) +4) _ Qo,i(uf A) ) Qi.r^+4) 



(^ M) +2) Qo,i(^ M) +2) Qr,!^) 



Qi,i («J 



(o,i) 2 



(8.29) 



8.2.4 Bethe equations for the Kac-Dynkin diagram <^) 1 

Using eq.(4.16) or eq.(4.15), we write the operator generating series: 
T(a,l,u + a+l) 2adu _ y ~_ lf 



Qho(u) e 2 du \ ( Qi,o(m) QiAu + 2) _ &w \ 1 / Q 2 .i{u) Qi,i(m + 2) _ &2d% 



(8.30) 



.Qi,o(« + 2) ;VQi,o( w + 2 ) Qi,i(*0 / \Q2,i{u + 2) Qi,i(u) 

In particular, 

T(l,l,u + 2) Qi,o(« + 4)0(u + 2) , Q ll0 (« + 4) Qi,i(«) (/>(u + 2) Qi,i(u) 



0(u + 3) Qi,o(« + 2)^(u + 4) Q 1 ,o(« + 2)Qi,i(u + 2)^(u + 4) Qi,i(u + 2)' 

Canceling the poles in eq.(8.31), we write the following Bethe equations: 

,(1,0)' 

1 



(8.31) 



h,l [Uj 



hi(u^-2 



«' 1) +2) Q 1 ,o(uf A) +2) 



(8.32) 



8.3 Examples of the integration algorithm 

To illustrate the general algorithm of section 6, we apply it to gl{2), gl{2\\) and gl{2\2) (super)algebras. 

8.3.1 gl(2) algebra 

In this case equation (6.36) reads 

r 2 ,o(l, s, u) =H-} Q (0,s, u) Ti j0 (l,a, u), < s < oo. (8.33) 

Substituting explicit expressions for H~2~ Q (0, s, u) (6.19) and T\$(\,s,u) (6.34), we immediately obtain the 
explicit formula for non-vanishing functions which do not belong to the boundaries: 

r 2>0 (l, s,u + s-l) = Q lt0 (u + 2s)Q 1 , (u - 2) £ - ^f^L+M. (8 . 34 ) 

^ QiM u + 2 J - 2)Qi,o{u + 2j) 
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8.3.2 gl(2\l) superalgebra 

In this case equation (6.36) reads 

T 2A (l,s,u) = H 2fi +(0,s,u) H~} Q (0,s,u) T lfi (l,s,u) = H 2 ,o+ (0, s, u) T 2fi (l,s,u) 7 < s < oo. (8.35) 

Substituting explicit expressions for H 20 + (0, s, u) (6.11) and T 2y0 (l,s,u) (8.34), we obtain, after straightfor- 
ward calculations, the following result for non-vanishing functions which do not belong to the boundaries: 

T 2 , 1 (l, Sl , +S -l) = Q 1 ,o(.+2 s) ( gMM + Ql , l( ,-2)g n . Q 2't"n + + o ■ + » 

\Qi,o{ u ) f^o Qi.o( u + 2j)Qi,o(« + 2j + 2) J 

The step function 0(s) = 1 at s > and otherwise. In the calculation, we have used eq. (3.45) to substitute 

Qi,o(u - 2)Q 2 ,i(") - Qi,o(«)Q2,i('U - 2) 



:,o(u) 



Qi,i(«-2) 



8.3.3 gl(2\2) superalgebra 

In this case equations (6.36) and (6.37) read 

T 2j2 (l,s,u) = H 2 ,i+(0,s,u) H 2 fi+(0,s,u) H~} Q (0,s,u) T lt0 (l,s,u) 

= H 2A +(0,s,u) H 2fi +(0,s,u) T 2<0 (l,s,u) = H 2 ^+(0,s,u) T 2 ,i(l,s,u), < s < oo (8.37) 

and 

T 2 , 2 (a, 1, u) = Hi+ 2 (a, 0, u) H + 2 (a, 0, u) H~ 2 _ (a, 0, u) T 0i i(a, 1, it), < a < oo . (8.38) 
Substituting explicit expressions for the shift operators, wc finally obtain: 

'Q 2 , 2 (u) , Qa,a(«)Qi,i(ti-2)Q 2) o(u + 2) 



T 2 , 2 (l, s, u + s - 1) = Qi, (u + 2s) I - — — + — — -— ^ — 1 —= — - — — — 9{s - 1) 

s-2 
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, 2 (u-2)Q a ,i(u + 2) _ , ,>^ Q2,o(w + 2j + 4) 



O^NM" + 2) " - " 1} + QlMu ~ 2) g Qi,o(« + 2T+2)Q;,o(« + 2 j - + 4) ^ " 2) ) (8 ' 39) 
and 

T 2j2 (a,l,u- a + 1 = Qoa 2a I - — — + — - — — — — — — — — 6 (a - 1 

v <2o,i( u ) Qi,2{ u )QoA u )Qo,i{ u - 2) 

+ _ g 4 , «(. - 2)) .(M0, 

Qi, 2 {u)Qo^(u - 2) ^ QoA u ~ 2j - 2)Q ,i(w - 2j - 4) / 

In the calculation, we have used eq. (3.45) to substitute 

Qi,i(u - 2)Q 2:2 (u) - Qi,i(u)Q 2j 2(u - 2) = Q li2 (« - 2)Q 2> i(«) 

and 

Qo,i(w)Qi,2(u + 2)-Q ,i(w + 2)Qi.2(u) ^ , . „* 
^ — ^ = Qi,i{u + 2). 

The functions T 2j2 (l,s,w) and T 2i2 (a, l,u) form a complete set of non-vanishing T-functions which do not 
lie on the boundaries of the fat hook domain for the case of the g/(2|2) superalgebra. 
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9 Discussion 



In this paper we have dealt with finite dimensional representations of gl(K\M) and the periodic spin chains 
based on the rational i?-matriccs. Wc believe that the method is powerful enough to incorporate various 
generalizations, such as extension to twisted and open spin chains, to infinite-dimensional representations of 
non-compact version of the symmetry group as well as to models with more general i?-matriccs, including 
various exotic ones, as the Hubbard i?-matrix or the recently constructed AdS/CFT S-matrix. We expect 
that the Hirota relation should be the same for all these problems, with only difference being in the boundary 
and analyticity conditions. 

For example, the extension to spin chains with twisted boundary conditions can be accomplished by 
replacing the polynomial Q-functions by "Bloch polynomials" , i.e., functions of the form Q(u) = Ae KU Ylj( u ~ 
Uj), where k is related to the twisted boundary condition. This leads to appropriate simple modifications in 
the Bethe equations. 

The extension to supcrsymmctric spin chains with trigonometric i?-matrices is also straightforward. The 
Hirota equation and the boundary conditions remain the same. The only change is again in the analytical 
properties of the T-functions: in the trigonometric case they are "trigonometric polynomials", i.e., finite 
products of sin(?7(w — Uj)) or sum (j](u — Uj)). A hypothetical generalization to the elliptic case is much more 
interesting. As far as we know, supersymmetric quantum spin chains with elliptic i?-matrices were never 
discussed in the literature. On the other hand, the bosonic case suggests [16] that the functional relations 
between commuting integrals of motion are given by the same Hirota equation. The solutions are sought in 
the form of "elliptic polynomials", i.e., finite products of Jacobi theta-functions 9(rj(u — Uj)\r). From this 
point of view, it would be interesting to analyze elliptic polynomial solutions to the Hirota equation with 
boundary conditions of the fat hook type. They might solve (as yet hypothetical) supcrsymmctric quantum 
integrablc models with elliptic i?-matrices. 

It is also important to elaborate, within the framework of Hirota equations, a more direct approach to 
spin chains in typical representations of superalgebras (the ones having a continuous label). One would like 
to understand them not only in the sense of analytic continuation with respect to the representation label 
(as we demonstrated in this paper by simple examples) but also to find their place on other levels of the 
construction. A characterization of solutions to the Hirota equation that are responsible for spin chains in 
typical representations would be of particular importance. 
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Appendix A. Superalgebras and their representations 

For completeness, we list here some basic objects and notation related to Lie superalgebras gl(K\M) and 
their representations. The main references are [19]- [22]. 
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We use the notation p{a) for the Z2-grading of the index a: p{a) = if a is bosonic and p(a) = 1 if a 
is fermionic. The same notation is used to denote the grading of any objects (vectors, operators,...) with 
definite parity. Objects with definite parity are called homogeneous. 

A matrix A with matrix elements Aj (we imply that they are usual numbers, i.e., p(Aj) = 0) is said to 
be even (odd) if p(i) + p(j) is even (odd) for all non-vanishing elements of A. The i?-matrix (2.1) is an even 
matrix. 

Let X, Y be two graded spaces, and let e^, fj be the corresponding (homogeneous) basis vectors. For 
any two vectors x = J2i x i e ii V = J2i Vifi we have 

x ® y = 5> lCl ) ® {jufi) = J2(-iy^x iVj {<* ® fj) , 

so the components of the vector x ® y in the basis ® fj are (— l) p ^' p ^ Xiy 3 -. The action of the operator 
A ® B in X <g> Y" is defined on homogeneous objects as 4 ® B (x <g) y) = (-l)^)^^^) g> B(y). Matrix 
elements of the tensor product of even matrices A*, Bp are: 

(A ® B)j^ = (-i)p( Q )(f( l )+p(J'))^. J B| . 

This rule explains the origin of the sign factors in the graded Yang-Baxter equation. 

The Lie superalgcbra gl(K\M) can be most transparently defined through its matrix realization: it is 
the set of block matrices 

«-(£{!) < A1 > 

such that A, B, C, D are respectively K x K, K x M, M x K and M x M matrices. The even subalgebra 
gl(K\M)o has B = C = 0, the odd subalgebra (?Z(i4T|M)i has A = D = 0. For homogeneous elements, the 
bracket is given by 

[g,9'\ = .9s'-(-i) p(9)p(9 Vs- 

Note that gl(K\M)o = gl(K) © gl(M). For elements g realized as above the supertrace is defined by 
stig = trA- tiD. 

A basis for gl(K\M) consists of matrices Eij with entry 1 at position and otherwise. A Cartan 
subalgebra of gl(K\M) is spanned by the elements En, i = 1, . . . , K + M. The set of generators of gl(K\M) 
consists of the En and the elements -E^i+i and Ei+ij, i = 1, . . . , K + M — 1. The space dual to the 
Cartan subalgebra is spanned by the linear forms : g h- > An (i = 1,...,K) and Si : g ^ Da (i 

1 M ), where 5 is given by (Al). Let us choose the basis in the space dual to the Cartan subalgebra to be 

ei, . . . , tfc, S\, . . . , Sm- On this space there is a bilinear form induced by the supertrace in the superalgebra: 

N<y) = % > = (#i|<y) = , = Sij . 

Even (bosonic) roots are e, — and Si — Sj (i ^ j), odd (fermionic) roots are ±(ej — (5j). There are 
several choices of simple root systems. The distinguished simple root system has the form ctj = — e.; + i for 
i = 1, . . . , K — 1, olk = £jf - <5i, ax+j = Sj — Sj+i for j = 1, . . . , M — 1. All other simple root systems are 
obtained from this one by reflections with respect to odd roots a with (a\a) = 0. 

Elements of the space dual to the Cartan subalgebra are called the weights. The weight is an expression 
of the form 

K M 

A = ^A i e i + ^A i 5 i . 
i=i j=i 

Let A be a Young diagram formed by a sequence of non-negative non-increasing integers: A = (Ai, A2, . . .), 
Ai > A 2 > . . . > 0. To such a diagram one assigns the weight A, = Aj (i = 1, . . . , K), Aj = Xj (j = 1, . . . , M), 
where Aj = max (A^ — K, 0) and A^ is the height of the j-th column of the diagram A. It is implied that 
Ajr+i < M. 

A Kac-Dynkin label is a sequence b\, b%, . . . , bf[, ■ ■ ■ , bn+M-i, where all bj except bx are non- negative 
integers while bx may be any real number. There is a one-to-one correspondence between finite-dimensional 
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irreducible representations (irreps) of the superalgebra gl(K\M) and the Kac-Dynkin labels [19]. We consider 
covariant tensor irreps of the superalgebra gl(K\M). One can assign the following Kac-Dynkin label to any 
highest weight A associated with a Young diagram A as above: 

hi = Xi — Ai+i , i = 1, . . . , K - 1 , 

b K = Xk + Ai , (A2) 

bj +K = Xj - Xj+i , i = 1, • • • , M - 1 . 

Therefore, one associates a tensor irrep of gl(K\M) to any Young diagram. (The diagrams containing a 
rectangular subdiagram with K + 1 rows and M + 1 columns correspond to vanishing representations, so 
such diagrams are illegal, similarly to diagrams containing K + 1 rows for gl(K).) However, for supcralgebras 
this correspondence is not one-to-one. Different Young diagrams may correspond to equivalent irreps (i.e., to 
the same Kac-Dynkin label). In particular, this is the case for rectangular diagrams 8 when M + n columns 
of K boxes are replaced by M columns of K + n boxes [21]. 

There is a large class of finite-dimensional irreps of gl(K\M) which cannot be associated with a Young 
diagram. Given an irreducible tensor representation with the highest weight A, there is a one-parametric 
family of finite-dimensional irreps with the highest weight A(c) = A + c^\. =1 ej, where c is a real parameter. 
This yields the Kac-Dynkin label with a non-integer bx = Xk + Ai + c. 

One distinguishes typical and atypical irreps [19] (in physical terminology, long and short irreps, respec- 
tively). An irrep with the highest weight A is atypical iff there exists at least one pair i = 1, . . . , K, 
j = 1, . . . , M, such that 

(A + p,e i -S j ) = 0, (A3) 

where 

K M 

p=- J2(K - M -2i + l)e t + - J2(K + M — 2j + 1)5, . 

t=l " 3=1 

For typical irreps, there is no such pair that (A3) holds. All irreps with 6jf > M are typical. All 

rectangular irreps corresponding to diagrams having a rows and s columns with a < K or s < M are 
atypical. 

Appendix B. Auxiliary linear problems for the Hirota equation 

Let r = t(pi,P2,P3) be a function of three variables. For brevity, we denote 

r(pi + l,p 2 ,p 3 ) := Ti, r(pi,p2 + 1,P3) ■= r 2 , r{pi + l,p 2 + 1,253) := r 12 , etc. 

Let a/37 be any cyclic permutation of 123. Consider the following system of three linear equations for a 
function tp = tp(pi,P2,P3)- 

e d - +A 7 ^^ ) V = e 9 "V, a/3 7 = 123, 231, 312, (Bl) 

T a Tp J 

where A Q are parameters and d a = d/d Pa . Their compatibility implies the Hirota equation for r: 

A1T1T23 + A 2 t 2 ti3 + A3T3T12 = . (B2) 

To see this, consider the first and the second linear equations. Their compatibility means that the 
difference operators 



-e,( e e,_ X3 ZnA and e^fe^ + A!^ 

V T lT2 / V T 2 T 3 ) 



8 For brevity, we call irreps corresponding to the Young diagrams of rectangular shape rectangular irreps. 
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commute. A simple calculation shows that their commutativity is equivalent to the condition that the 
function 

A1T1T23 + A3T3T12 . . 

:= h(pi,p 3 ) 

T2T13 

does not depend on p2- Compatibility with the third linear problem implies that the function h must be 
a constant equal to —A3, whence the Hirota equation follows. In the case when the function h can be 
hxed in some other way (for example, from the boundary conditions), just two linear problems arc enough 
to represent the Hirota equation as a discrete zero curvature condition. Moreover, the specific form of the 
coefficient functions in the difference operators (Bl) implies that the compatibility follows from the existence 
of just one non-trivial common solution (cf. [44]). 

In terms of the function tp — tpr the linear problems acquire the form 

T-ytpp — Tpip 7 + \ a Tp 1 ip = , a/37 = 123, 231, 312 . (B3) 
From the first and the second ones we have 

7"3^2 + A1T23V? T 3 <y9l - A 2 T13<£ 
T 2 = , Tl = . 

Plugging this into the Hirota equation, we obtain another linear problem compatible with the previous ones: 



AlT23<£>l + A2Ti3< ( 92 + A 3 Ti2<^3 = . 

The four linear problems can be combined into a single matrix equation as follows [43] : 



(B4) 



/ 





T3 


-72 


A1T23 


\ 


/ 


<P1 




-T3 





Tl 


A 2 Ti3 






p2 




To 







A3T12 






¥>3 


V 


-AlT 23 


— A 2 Tl3 


— A3T12 





/ 


V 


<P 



(B5) 



The determinant of the antisymmetric matrix in the left hand side is equal to (A1T1T23 + A2T2T13 + A3T3T12) 2 . 
It vanishes if r obeys the Hirota equation, and the rank of the matrix is 2 in this case, so only two of the 
four equations are linearly independent. 

One may regard system (B3) as linear equations for r with coefficients p. Shifting the variables pp — > 
pp — 1, p 7 — ► pj — 1, and then passing to the new variables pi.2,3 — * — Pi, 2,3, one sees that the form of this 
system remains the same. Since the Hirota equation is invariant under the simultaneous change of the signs 
of all variables, the compatibility condition is the same Hirota equation for tp: 

Al<^23yi + A2<£l3<£2 + A3V12V3 = 0. (B6) 
Let us pass to the "laboratory" variables xi,X2,a;3 according to the formulas 

Pi = \{-EiXi+S2X2+£3Xz) 
p 2 = \{£lXi - £2X2 + £3X3) 
P3 = ^(£1X1 + £2X2 ~ £3X3) , 

where £ a = ±1 is a fixed set of signs (clearly, there are 2 3 = 8 possible choices). The inverse transformation 
reads 

xi = £i(p2 +P3) , x 2 = ez(pi +P3) , x 3 = £ 3 (pi +p 2 ) ■ 
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We introduce the functions T(x\, x 2 , X3) = r(pi,p2,P3) and F(xi,X2, X3) = f{p\,P2,P3), where the variables 
x a and p a are connected by the formulas given above. In the "laboratory" variables, the system of four 
linear problems takes the form 



/ 





T12 


-T13 


A1T1123 


\ 


/ F23 


\ 




-T12 





?23 


A2? 1 i223 




F13 






T13 


— T23 





A3T1233 








V 


— A1T1123 


— A2T1223 


— A3T1233 





/ 


V 


/ 



= 0. 



(B7) 



where we denote T\ = T(xi+ei,X2,x 3 ), T 12 = T{x\+e\, x 2 +s 2 , x 3 ), T n2 3 — T(x 1 + 2ei, x 2 +£2, x 3 +£ 3 ), etc 
(and similarly for F). 

The compatibility of these linear problems implies the Hirota equation 

A1T1123T23 + A2T1223T13 + A3T1233T12 = . 
Shifting the variables (x a — > x a — e a ), we get the equation 

XlT(xr\€i, X 2 , X3)T(Xi-ei,X2,X3)+\2T(xi,X2+£2,X3)T(xi,X2-e2, X3)+\3T(xi,X2,X 3 i€3)T(x 1 ,X2,X3-£3) = 

Note that it is the same for any choice of the e's. Linear equations (B7) provide inequivalent Backlund 
transformations for it. However, only four of them (corresponding, say, to the choices E\ = £2 = £3 = 1, 
— £1 = £2 = £3 = 1, £1 = —£2 = £3 = lj £1 = £2 = —£3 = 1) are actually different because the simultaneous 
change of signs of all e's means passing to the "conjugate" system of linear problems, where the roles of T 
and F are interchanged. Equation (3.13) in the main text corresponds to the choice x\ = a, X2 = s, X3 = u, 
Ai = A2 = —A3 = 1, £1 = —£2 = £3 = 1. 



Appendix C. On bilinear equations (3.38)- (3.43) and their compat- 
ibility with (3.29), (3.30) 

Here we complete the proof of equation (3.38) and present some more details on compatibility of bilinear 
equations (3.29), (3.30) and (3.38)- (3.43). 

As we have seen, the pair of equations (3.33), (3.34) implies the relation 

Tk, m (a, 8+1, u)T k+ x >m+ i{a, s, u+1) - T k>m (a, s, u+l)T fc+ i )m+ i(a, 8+1, u) 
= fk,m(a, u+s) T k+X>m {a, s, u+l)T k>m+ i(a, s+1, u) , 



(CI) 



where fk.m(o-,u + s) is an arbitrary function of k,m and a, u + s. In the same way, the pair of equations 
(3.35), (3.36) implies the relation 



(C2) 



T ktm (a-1, s, u)T k+ x >m+1 (a, s, u + 1) - T k>m (a, s, u + l)T k+ i^ m+ i(a-l, s, u) 
= 9k,m{s,u-a)T k>m+1 (a-l, s,u)T k+ i >m (a, s,u+l) , 
where g k>m (s, u~a) is an arbitrary function of k, m and s, u — a. On the other hand, we have the identity 
T k ^ n (a,s + l,it)T fc+lim+ i(o, s,u+l)- T k>m (a, s,u+ l)T k+ i. rn+1 (a, s + l,u) 



Tfc, m+ i(a, s + 1, u)T k+1:m (a, s, u + 1) 

Tk, m (a ^ 1, s,u)T k+1 , m+1 (a, s,u+l) - T fc;m (a, s,u+ l)T k+1 , m+1 (a - l,s,u) 
T kym +i(a - l,s,u)T k+li7n (a, s,u+ 1) 



(C3) 
(C4) 
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which is straightforwardly proved by passing to the common denominator, grouping together similar terms 
and using equations (3.29), (3.30). Its left hand side is /fc, m (a,u+s), and the right hand side is gk,m(St u— a). 
Therefore, we conclude that 

fk, m (a, u + s) = g k , m (s, u-a) . 

This equality is possible only if fk.m(o-,u + s) actually depends on the difference of its arguments, i.e., on 
u + s — a. This fact has been used in section 3.4 to prove that /fe.,„(a, u + s) = 1. 

The system of bilinear equations (3.38)-(3.43) can be represented in the form 



(e 



<>u 



(e 



d a 



r e n) 
-e 9 °) 





S U) 


- 1 k,m+l V u 


i 4- 


1 




Tk,m{a, s, 


u) 
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1, 


u) T k , m {a, S,U + 1) 




s, u) 




s, u 




1) T k+ i, m (a + l,s,u) 


Tk, m (a, s, 


u) 




s, u 




1) T k ,m(a + l,s,u) 


Tk+l, m +l( a 


s, u) 




s + 


1. 


u) T k+liTn (a + l,s,u) 


Tk, m (a, s, 


u) 




s + 


1. 


u) T kt m(a + l,s,u) 


Tk+\,m+i{a 


s, u) 


2\m+i(a 


- 1, 


s 


+ l,u- 1) T fc+ i !m (a + l,s,u) 


Tk, m (a, s, 


u) 




n(a, 


s, 


u- 1) Tfe jm (a, s+ l,it) 




s, u) 




s + 




it) T fc+ i ;m (a + l,s- l,u + 1) 


Tk, m {a, s, 


u) 


Tk, 


n(a, 


s, 


it + 1) T fc . m (a + l,s,u) 


7fc+l,m+l( a 


s, u) 




-I, 


s 


+ l,u- 1) T fc+lim (a,s,u+ 1) 




u) 








l,s,u) r fei7n (a, s + l,u) 



(C5) 
(C6) 
(C7) 
(C8) 
(C9) 
(CIO) 



These equations have a similar structure: different commuting difference operators in the left hand sides 
act on the same fraction of the T- functions. It is a simple exercise to verify that these equations are self- 
consistent if equations (3.29), (3.30) are satisfied. The self-consistency conditions obviously follow from 
the commutativity of the difference operators in the left hand sides. Now, let us demonstrate that all the 
equations from (3.29), (3.30) are encoded in the system (C5)-(C10). To this end, we subtract the sum of 
equations (C5) and (C9) (equations (C7) and (C8)) from eq. (C7) (eq. (C6)) and find that the resulting 
equation reproduces the first equation from (3.30) (respectively, the second equation from (3.29)). The 
remaining second (first) equation from BT2 (3.30) (BT1 (3.29)) results from the consistency condition of 
equations (C6) and (C8) (respectively, equations (C5) and (C9)). 

The connections between the different bilinear relations for the transfer matrices Tfc im (a, s,w) become 
more transparent if one represents them in a matrix form. Namely, two pairs of equations (3.38), (3.42) and 
(3.39), (3.41) can be rewritten as 



Tk, m (a, s+1, u) 
-Tfc, m (a+1, s, u) 

Tk, m +i(a,s + l,u) 



-7\ + i, m+ i(a, s + 1, it) 

Tfcfi >m+ i(a+l, s, u) 
Tk+i, m (a,s,u+l) 

Tfe+i, m (a + l,s-l,u+l) 



7fc+i, m+ i(a, s, u+1) 
Tk, m {a, s, u + 1) 



(Cll) 



and 



-Tk, m (a, s+1, u) T fe+ i !m+ i(a, s+1, u) 

-T k: m(a+l, s, u) T k+1:m+ i(a+l, s, it) 

T k)7rH1 (a-l, s+l,it-l) 

7fc, m+ i(a, s, u-1) 



Tk+i, m +i(a, s, u-1) 
T kt m(a,s,u-1) 



T fc+ i !ln (a+l,s,u) 



(C12) 
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respectively. Then, multiplying both sides of equations (Cll) and (C12) by the matrices inverse to the ones 
in their left hand sides (for the calculation of the determinant eq. (3.40) is used), the resulting equations 
reproduce respectively the first and the second equations from the linear systems BT1 and BT2 (3.29), 
(3.30). Similarly, representing the two equations of the last line of the matrix equation (3.32) in the matrix 
form 

-7fc+i,m+i(a-l,s,u) T k+ i^ m+ i(a, s, u—1) \ / T fc+1 , m (a + l, s, u) 
-T ktm (a-l,s,u) T k , rn (a,s,u-l) ) \ T k+1 j n (a, s, u+1) 

(T fc+ i. m+ i(a, S+l,u) \ 
(C13) 

multiplying the both sides by the matrix inverse to the one in the l.h.s., and using eq. (3.42) to calculate the 
determinant of this matrix, one arrives at equations (3.41), (3.43). 



Appendix D. An alternative derivation of the Bethe equations 

Here we give a direct derivation of the Bethe equations from the pair of linear problems (Backlund transfor- 
mation) (3.29). It does not use the Hirota equation for the Q-functions. 

We know that we can use the k (respectively, m) flow of the BT1 (resp., BT2) to decrease rank K (resp., 
M) of the original problem to K — 1 (resp., M — 1). For the step by step passing from (K, M) to (0, 0) we can 
choose any zigzag path of the type drown in Fig. 17. The interior boundaries of the domain of non-vanishing 
T's move towards the exterior ones until the domain collapses to the horizontal and vertical lines, so that 
the original problem gets "undressed". The final solution can be formulated in terms of Bethe equations. 
Let us derive them for the simplest path with just one turn. 

First, using the transformation BT1 at each step, we move the horizontal interior boundary a = K, s > M 
to the half-line a — 0, s > M, as shown in Fig. 10. Consider the second equation from eq.(3.29) at s = 
(Fig. 12, position 4): 

Tk,M(a, ^,u+l)Qk-i,M{u+a) - Qk,M(u+a)T k _ hM (a, l,u+l) = Qk,M(u+a +2)T fe _ M/ (a-l, l,u), (Dl) 
where k = 1, . . . , K , or, denoting P k (u) = T fejM (a, l,u - a + 1), Qk,u(u) = Qk{u): 

PS(u)Q k -i(u) - P c U{u)Qk{u) = P^Zl{u - 2)Q k (u + 2), k = 1, . • • , K. (D2) 
We know that P k {u) and Q k {u) should be polynomials in u. We set: 

Q*(tt)=n(«-«? ) ) . (D3) 

Taking equation (D2) at zeroes of Qk-i(u), Q k (u) and Q k (u + 2), we obtain the equations 

( u[ k - r) + 2 



p" 



(««- 2) -2) = PU(nf-2)Q k (^-2). (D4) 



Dividing the second equation (with the shift a — > a + 1) by the third one and excluding the ratio of P's with 
the help of the first equation, we arrive at the following standard system of Bethe equations: 



fcfi (u?>+2) Q t (uf-2) Q t -r(uf) 
?t + i \uf) Q t («f+2) Qt-i (uf-2) 



-1 (j = l,...J t , t=l,...,K-l). (D5) 
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Then, using the transformation BT2 eq.(3.30) at each step, we move (as shown in Fig. 14) the vertical 
interior boundary s = M, a > to the half-line s = 0, a > 0. Consider the second equation of (3.30) 
at s = (Fig. 12, position 4). In a similar way, we get the Bethc equations for the roots Uj~ M ^ of 

Qo,m(u) := Qm-At(u): 



(««-2) 0t _ 1 („« + : 



t+ i [uf -2)Q t [uf +2 \Q t _ x 



-1 (j = l,...J t , < = -!,..., -M+l). (D6) 



To get the missing equation for zeros Uj of the function Qo{u) — Qo.m( u )- we use the second equation of 
BT1 at k = 1, m = M, and the second equation of BT2 at A: = 0, m = M: 

P?(u)Q (u)-Pg (u)Qi(u) = P °- 1 (u-2)Q 1 (u + 2), 

iZ?(u)Qo(«) - P tt («)Q-i(«) = P a_1 (^-2)g-i(w + 2), (D7) 

where -Pf (u) = Ti m (Oj l,u — a + 1), Rf(u) = ?b Af— 1, u — a + 1). At the "fermionic" roots of Qo( u ) = 



J the ratio of these two equations gives the missing system of Bethe equation: 



(0) 



Qi( 



«f +2 Q-r M« 



7 I A;.- (D8) 



The "boundary conditions" for the system of nested Bcthe equations are: Qk{u) = Qk,m{u) = <fi( u ) ( a 
given polynomial function), Q-m(u) = Qo,o(u) = 1. 

Note that our numbering of the Q's here differs slightly from the one in section 4.5: the Q's here are 
numbered from — M to K while the same Q's there have the numbers from to K + M. Taking this into 
account, it is easy to see that the chain of Bethe equations obtained here coincides with (5.8), where the 
signs are chosen as p 1 = p 2 = . . . = pm = -1, Pm+i = Pm+2 = ■ ■■ = pm+k = 1. 

This is only one of possible sets of Bethe equations for the same model. We could also apply the 
elementary moves BT1 and BT2 in any other order. At each change of the direction we obtain an equation 
of the fermionic type (D8). The different sets of equations lead to the same solution for the T- functions 
T(a, s, u) = Tk,m{cl, s, u) at the highest level of the hierarchy. The different systems of Bcthe equations are 
known to be related by certain "duality transformations" [14, 45]. As we have seen in the main text, they 
admit a natural explanation in terms of a discrete "zero curvature condition" on the (k,m) lattice. 

Appendix E. Comparison of the results for gl(2\l) algebra with the 
results of [6] 



Let us compare eq.(6.1) of [6], with the notation fi = — ^u, = \u^\ = , b = s + : 

Ju Jv n 

Qi(u) = A ( u - u ?) ' Q2{u) = \{{u-uf), 0±(«) = JJ(«±r i ), 
j'=i j=i i =1 
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with our eq. (8.23). The Baxter TQ-relation from [6] for the (b, 1/2) irrep in the auxiliary space and 1/2) 
irreps at the sites of the chain now reads as follows: 



r s|{r ( }0) 



Q 2 (u + s) 
Q 2 {u - s) 

i(u - s + 1) Q 2 (u + s) <t>+{u 



Qi(u — s — 1) Q 2 (u — s) <f)-(u — s) 

Qi(u — s — 3) Q 2 (u + s) 4> + (u-s) 
Qi(u — s — 1) Q 2 (u — s — 2) <t>-{u - s) 

Qiju + s) - s) (j>+(u -s-2) 

Q 2 (u-s-2) 4>-{u-s) <j)-(u-s-2) ' { ' 



We see that it coincides with our eq.(8.23) up to the redefinition of (^-functions. 
Let us give also the Bethe equations ensuring the polynomiality of T s \{ ri y(u): 



{2 

1. Canceling poles at u = u K 



j 1 ' 

( 2 h n.fJ 2 ) 



(E2) 



<t>-{uy>) _ Qi(uy> + i) 

Muf) Qi(uf-i) 

we arrive at the Bethe equations for the fermionic node. 

2. Canceling the poles at u = lip + s + 1, 

i _ Q 1 {uf+2)Q 2 (uf-l) 

we get the Bethe equations for the bosonic node. (We note a mistake in the denominator of the r.h.s. 
of eq. (6.3) of [6].) 

(2) 

3. The condition of canceling poles at u = Uj ' 
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